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Abstract. The cellular automata (CA) approach to traffic modeling is extended 
to allow for spatially homogeneous steady state solutions that cover a two 
dimensional region in the flow-density plane. Hence these models fulfill a basic 
postulate of a three-phase traffic theory proposed by Kerner. This is achieved by a 
synchronization distance, within which a vehicle always tries to adjust its speed to 
the one of the vehicle in front. In the CA models presented, the modelling of the 
free and safe speeds, the slow-to-start rules as well as some contributions to noise 
are based on the ideas of the Nagel-Schreckenberg type modelling. It is shown 
that the proposed CA models can be very transparent and still reproduce the two 
main types of congested patterns (the general pattern and the synchronized flow 
pattern) as well as their dependence on the flows near an on-ramp, in qualitative 
agreement with the recently developed continuum version of the three-phase traffic 
theory [B. S. Kerner and S. L. Klenov. 2002. J. Phys. A: Math. Gen. 35, 
L31]. These features are qualitatively different than in previously considered CA 
traffic models. The probability of the breakdown phenomenon (i.e., of the phase 
transition from free flow to synchronized flow) as function of the flow rate to the 
on-ramp and of the flow rate on the road upstream of the on-ramp is investigated. 
The capacity drops at the on-ramp which occur due to the formation of different 
congested patterns are calculated. 
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1. Introduction 

Traffic on a highway can be either free or congested. In empirical investigations 
congested traffic shows a very complex spatial-temporal behaviour (see the reviews [1, 

2, 3, 4]). Based on a recent empirical study [5, 6] Kerner found out that in congested 
traffic two different traffic phases should be distinguished: "synchronized flow" and 
"wide moving jam". Therefore, there are three traffic phases: 1. free flow, 2. 
synchronized flow, 3. wide moving jam. 

A wide moving jam is a localized structure moving upstream and limited by two 
fronts where the vehicle speed changes sharply, i.e. within a region that is small 
compared to the distance between the fronts. A wide moving jam propagates through 
either free or synchronized flows and through any bottlenecks (e.g. at on-ramps) 
keeping the velocity of its downstream front [4, 7]. In this respect it differs from 
synchronized flow, the downstream front of which is usually fixed at the bottleneck, 
where it occurred. Such empirical spatial-temporal features of "wide moving jams" and 
"synchronized flow" are the basis for the distinction of these traffic phases in congested 
traffic rather than a behaviour of traffic data in the flow-density plane [4, 7, 8]. 

Wide moving jams do not emerge spontaneously in free flow (with the possible 
exception that synchronized flow is somehow prohibited [9]). Instead, there is a 
sequence of two first order phase transitions [6] : First the transition from free flow to 
synchronized flow occurs (called the F-^S-transition), and only later and usually at a 
different location moving jams emerge in the synchronized flow (the S— > J-transition) . 

Different explanations of these empirical findings have been proposed by various 
groups in the last years, but so far they remain controversial (see e.g., [10, 11, 12, 13, 
14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26] and the excellent review by Helbing [1]). 

1.1. The fundamental diagram approach and cellular automata models of the 
Nagel-Schreckenberg (NaSch) type 

Empirical observations show that the higher the vehicle density is the lower the average 
vehicle speed. The average flow rate, which is the product of the average speed and 
the density, is a function of the density which has a maximum. This curve in the 
flow-density plane is called the fundamental diagram [1, 2, 27]. 

Apparently the empirical fundamental diagram was the reason that already the 
first traffic flow models [28, 29, 30, 31, 32] were based upon the postulate that 
hypothetical spatially homogeneous and time-independent solutions, where all vehicles 
have the same distances to their neighbours and move with the same constant speed, 
(steady state solutions for short) exist that are related to a fundamental diagram, 
i.e. a curve in the flow-density plane. (Steady state model solutions are often called 
"equilibrium" solutions or "equilibrium" states of a traffic flow model. In this paper we 
will use the term "steady states".) A subset of these steady states would be unstable 
with respect to noise or external perturbations. This postulate underlies almost all 
traffic flow modeling approaches up to now [1, 2] in the sense that the models are 
constructed such that in the unperturbed, noiseless limit they have a fundamental 
diagram of steady states, i.e. the steady states form a curve in the flow-density plane. 

The congested patterns which arc calculated from these models for a homogeneous 
road (i.e. a road without bottlenecks) [1, 33, 34, 35] are due to the instability of steady 
states of the fundamental diagram within some range of vehicle densities. When 
perturbed, they decay into a single or a sequence of wide moving jams ("moving 
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clusters" [33]), whose outflows accelerate to rather high speeds comparable to free flow. 
This is why we find it helpful to classify these models as belonging to what we call the 
"fundamental diagram approach" . In the next subsection another class of models will 
be described, which belong to the "three-phase traffic theory" and lead to qualitatively 
different congested patterns, which are in better agreement with real observations. Its 
hallmark is a third kind of flow in addition to the wide moving jams and free flow, 
in which all vehicles interact strongly but nevertheless drive smoothly at a reduced, 
rather uniform speed (synchronized flow). In view of the rich phenomenologics of the 
two model classes in particular near an on-ramp, the distinction of the two classes is 
admittedly a bit simplified at this point, but a more detailed comparison will be given 
in Sec. 5. 

Apart from the congested patterns on a homogeneous road the dynamical 
behaviour near an on-ramp differs significantly in the fundamental diagram 
approach [1, 10, 11, 12, 13, 14, 34, 15] and in the three-phase traffic theory [8, 24]. 
It depends on both the flow rate to the on-ramp q on and the initial flow rate on the 
road upstream of the on-ramp, <?;„. Therefore it can be conveniently characterized by 
a so-called diagram of congested patterns. This is a map with coordinates q on and q m 
of the regions, where different congested patterns upstream of the on-ramp occur. 

In this paper we focus on cellular automata (CA) traffic flow modeling which 
was pioneered by Nagel and Schreckenberg in 1991/1992 [36]. The original Nagcl- 
Schreckenberg model (NaSch model) as well as all subsequent modifications and 
refinements of it [22, 23, 36, 37, 38, 39] belong to the fundamental diagram approach 
(Fig. 1 (a)). See also the reviews by Wolf [40], Chowdhury et al. [2] and Hclbing [1]. 
In this paper CA models belonging to the three-phase traffic theory will be proposed, 
and it will be shown that their congested patterns differ qualitatively from the ones 
in the CA models of the NaSch-type. 

Within the fundamental diagram approach there is a subset of traffic models 
which show spatially homogeneous high density states with low speed. At first sight 
these states look like "synchronized flow" , but this interpretation is incompatible with 
the observed dynamical behaviour near an on-ramp, as will be discussed now. In the 
diagram of congested patterns near an on-ramp obtained by Helbing et al. [1, 10] such 
high density states with low speed occur upstream of the on-ramp, if the flow rate 
q on to the on-ramp is high enough. Hclbing et al. called these states "homogeneous 
congested traffic" (HCT) [10] and proposed to identify it with synchronized flow [1]. 
In HCT no moving jams occur spontaneously. 

By contrast, in empirical observations of synchronized flow at low vehicle speed 
and high enough vehicle density moving jams do emerge spontaneously. These moving 
jams are particularly likely to occur, when due to the high flow rate to the on-ramp 
a strong compression of synchronized flow appears ("the pinch effect") [6, 8]. In fact, 
moving jams are only observed to emerge spontaneously in synchronized flow upstream 
of an on-ramp on which the flow rate is high enough [8]. At lower flow rates to the 
on-ramp synchronized flow of higher vehicle speed can occur in which the nuclcation 
of moving jams was not observed. 

1.2. The three-phase traffic theory 

To explain the above and other empirical results, Kerner introduced a three-phase 
traffic theory which postulates that the steady states of synchronized flow cover a 2D 
region in the flow-density plane, i.e. there is no fundamental diagram of traffic flow 
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in this theory [3, 6, 41, 42, 43]. 

This is not excluded by the empirical fact mentioned above, that a given vehicle 
density determines the average vehicle speed. Even if there is a continuous interval of 
different vehicle speeds at the same distance between vehicles (at the same density) , as 
is indicated by empirical observations of synchronized flow, obviously their averaging 
leads to one value at the given density. A 2D region of steady states in the flow-density 
plane is also not ruled out by car following experiments, where a driver has the task 
to follow a specific leading car and not loose contact with it (e.g., [44]). In such a 
situation, the gap between the cars will be biased towards the security gap depending 
on the speed of the leading car. In synchronized flow the situation is different: The 
gap between cars can be much larger than the security gap. 

As it has recently been postulated on very general grounds [8, 45] and 
demonstrated for a microscopic traffic model [24] , the fact whether the steady states 
of a mathematical description of traffic belong to a curve or to a 2D-region in the 
flow-density plane qualitatively changes the basic non-linear spatial-temporal features 
of the congested patterns which the model allows. 

In particular the diagram of the congested patterns at on-ramps is qualitatively 
different in three-phase traffic theory [8, 45, 24] from the diagram obtained in the 
fundamental diagram approach [1, 10, 12, 13, 14, 15, 17]. Specifically, at a high 
flow rate to the on-ramp, instead of HCT without moving jams, we find that moving 
jams always spontaneously emerge in synchronized flow of low vehicle speed and high 
density, whereas synchronized flow of higher vehicle speed can exist for a long time 
without an occurrence of moving jams [8, 24, 45]. This agrees with the empirical 
observations [8]. 

The microscopic model proposed by Kerner and Klcnov in order to derive the 
diagram of congested patterns within the three-phase traffic theory is relatively 
complex [24] . The main aim of this article is a derivation of cellular automata models 
within the three-phase traffic theory (where steady states of the models cover 2D- 
regions in the flow-density plane) which on the one hand are the most simple ones 
and on the other hand are able to reproduce the diagram of congested patterns found 
in [8, 24, 45]. The article is organized as follows. First, in the next section several 
cellular automata models with qualitatively different 2D-regions in the flow-density 
plane for steady states will be introduced. Second, the congested patterns at an on- 
ramp and their diagrams will be numerically derived for these models and compared 
with one another and with the results in [24]. Third, the probability of the breakdown 
phenomenon (the F^S transition) at an on-ramp is studied and the capacity drop is 
calculated. Finally, the results of CA-models in the three-phase traffic theory and in 
the fundamental diagram approach are compared. 

2. Cellular automata models within three-phase traffic theory 

2.1. Equations of motion 

The starting point of CA-modeling of three-phase traffic theory is the basic set of 
rules from [24] which provides a 2D-region of steady states. Denoting the speed and 
space coordinate of a vehicle at discrete time t = nr, n = 0, 1,2, .. by v n and x n , 
respectively, the basic rules are rewritten for CA-models in the form: 

v n+ i = max(0,min(v frco , v Si „, v Ci „)), x n+1 = x n + v n+1 r. (1) 
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Wf roo is the maximum speed of the vehicles. It is assumed to be the same for all vehicles 
in this paper. w s „ is the save speed which must not be exceeded in order to avoid 
collisions. In general it depends on the space gap between vehicles, g n — x^ n — x n — d 
and the speed ve, n [46, 40], where the lower index I marks functions (or values) related 
to the vehicle in front of the one at x n , the "leading vehicle" , and al is the vehicle length 
(assumed to be the same for all vehicles in this paper). For the sake of comparability 
we neglect the W£ jn -dependence and choose the same expression as in the standard 
NaSch-model [36]: 



The crucial difference compared to previous CA-models is that the acceleration 
behaviour given by v Cj „ (the rule of "speed change") depends on, whether the leading 
car is within a "synchronization distance" D n = D{v n ) or further away [24]. At 
sufficiently large distances from the leading vehicle, one simply accelerates with an 
acceleration a, which is assumed to be the same for all vehicles and independent of 
time in this paper. However within the synchronization distance the vehicle tends to 
adjust its speed to the one of the leading vehicle, i.e. it decelerates with deceleration b 
if it is faster, and accelerates with a, if it becomes slower than the leading vehicle. The 
deceleration b should not be confused with braking for safety purposes (i.e. in order 
not to exceed v s ^ n ). In practice the speed adjustment within the synchronization 
distance can often be achieved without braking at all simply as a result of rolling 
friction of the wheels with the road. Therefore the deceleration b is typically smaller 
than the braking capability of a vehicle. For simplicity we set b = a in this paper, so 
that the speed change per time step within the synchronization distance is given by 
Av n — ar sgn(u^„ — v n ) , where sgn(x) is 1 for x > 0, for x — and —1 for x < 0. 
In summary, 



We want to emphasise that this rule decouples speed and gap between vehicles for 
dense traffic. This can be seen by assuming that vehicles drive behind each other 
with the same speed v. According to (3) neither the speed nor the gaps will change, 
provided, all the distances are anywhere between the safe distance d + vt and the 
synchronization distance D(v) > d + vr. There is neither a speed-dependent distance, 
which individual drivers prefer, nor is there a distance-dependent optimal speed. 
This is the principal conceptual difference between three-phase traffic theory and the 
fundamental diagram approach, and it is the reason, why in three phase traffic theory 
the steady states fill a two dimensional region in the flow-density-plane, while in the 
fundamental diagram approach they lie on a curve. 

Let us contrast (3) with two models belonging to the fundamental diagram 
approach, the NaSch-model with "comfortable driving" recently put forward by 
Knospe et al. [22, 23, 47] on the one hand, and Wiedemann's modelling approach 
[48] on the other. Knospe et al. [22, 23, 47] put forward an extension of the NaSch 
CA-modcl, in which a driver starts to brake within some interaction horizon, as soon 
as he sees the brake lights of the vehicle in front being switched on. If nobody puts 
on the brakes, vehicles would close up to the safe distance, which is a function of the 
speed. In this sense, the term "comfortable driving" is a bit misleading: The behaviour 
modeled is more accurately described as "comfortable braking" . In contrast, the speed 
synchronisation (3) in our model happens always, whether someone brakes or not. It 



v s , n = g n /r. 



(2) 




for x iiTl - x n > D 
v n ) for x t , n -x n < D 



(3) 
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reflects what we believe to be the typical way, in which drivers take into account the 
vehicle in front of them. 

In Wiedemann's modelling approach [48] each vehicle has a prefered following 
distance, but convergence to it is hindered due to imperfect perception, so that the 
actual distance has an oscillatory behaviour. This model belongs to the fundamental 
diagram approach, too, as the steady state solutions have a unique relationship 
between speed and distance between vehicles. 

Of course, in models with a fundamental diagram of steady states, fluctuations 
and external perturbations let the system evolve in time through a 2D region in the 
flow-density plane as well. However, the dynamics is governed locally by steady state 
properties, the unstable steady states acting as "rcpcllors" and the stable ones as 
"attractors" . If the steady states form a 2D region, part of which is stable and part 
of which is metastable, as is the case in three-phase traffic theory, the dynamics is 
fundamentally different. This leads also to qualitative differences between the patterns 
of congested traffic obtained in three-phase traffic theory or in the fundamental 
diagram approach, respectively, as will be shown in detail below. 

2.2. Synchronization distance 

The conditions (1), (3) are the basis of the cellular automata models under 
consideration. It will be shown that this allows different formulations for fluctuations, 
acceleration, deceleration and for the synchronization distance D n which all lead to 
qualitatively the same features of congested patterns and the same diagram of these 
patterns as postulated in three-phase traffic theory [8, 45] and in agreement with the 
continuum model of Kcrncr and Klenov in [24] . 

In particular, let us consider two different formulations for the dependence of 
the synchronization distance D n on the vehicle speed. In the first formulation, the 
synchronization distance D n in (3) is a linear function of the vehicle speed: 



In the second formulation, the synchronization distance D n in (3) is a non-linear 
function of the vehicle speed: 



In (4) and (5) di, k and (3 are positive constants. Both formulations lead to 2D-regions 
of steady states in the flow-density plane. 

2.3. Steady states 

Whereas in models belonging to the fundamental diagram approach (e.g., [10, 12, 16, 
18, 22, 23, 36, 39, 46, 49, 47] and the reviews [2, 1]) a vehicle would close up to the 
leading one adjusting its speed and gap as required by secure driving, in models with 
the basic structure (1), (3), a driver within the synchronization distance D n adapts 
his speed to the one of the vehicle in front without caring, what the precise gap is, as 
long as it is safe. This explains why there is no unique flow-density relationship for 
steady states in the present CA models. 

In steady states all accelerations must be zero. Then the time-index n can be 
dropped in the above formulas. According to (1) - (3), there are two possibilities: 
Either the synchronization distance D{v) < g + d and the speed is v = Vf ICC} or 



D{v n ) = dx + kv n r. 



(4) 



D{v n ) = d + v n r + j3vl/2a. 



(5) 



D(v)>g + d and v = v e < min(w frco , v s (g, v j). 



(6) 
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Thus, in steady states all speeds are equal v. The conditions v has to fulfill are only 
equal, if also the gaps g are all equal. Therefore we defined steady states above as 
time-independent and homogeneous. 

The density p and the flow rate q are related to the gap g and the speed v by 

p=l/(x e -x) = l/(g + d), q = pv = v/(g + d). (7) 

Because v and g are integer in CA-models, the steady states do not form a continuum 
in the flow-density plane as they do in [24]. However, the inequalities of (6) define 
a two-dimensional region in the flow-density plane, in which steady states exist. As 
in [24] it is limited by three boundaries (Figs. 1 (b) and 2 (a, b)), the upper line U, 
the lower curve L, and the left line F. The parameters of the lines F and U are chosen 
to be the same for all CA models under consideration. Note that without the lower 
boundary L, the lines F and U constitute the fundamental diagram of the NaSch 
CA-model (Fig. 1 (a)) [36]. 

The left boundary F is given by q = pvf rcc . This is free flow, where the flow rate 
q is not restricted by safety-requirements. On the upper boundary U the flow rate is 
determined by the safe speed v s . For example, inserting (2) and (7) it is given by 

q = pv s = {l- pd)/r. (8) 

The lower boundary L is determined by the synchronization distance D: A steady 
state with density p and a speed v < v hcc requires that D(v) > 1/p with equality on 
the lower boundary L. For example, using (4) with d\ = d one obtains (see Fig. 1(b)) 

q = (1 - pd)/kr. (9) 

In the second model (5) the lower boundary L is a non-linear curve (see Fig. 2(a)): 

q= r{\l 1 + ~p {1 ~ pd) ~ ^ ~ P=E T' (10) 

This curve has the upper line U as a tangent at p = p ma x = 1/d. As will be shown 
below, this allows to reproduce qualitatively the diagram of congested patterns of 
three-phase traffic theory with simpler fluctuations than what is needed in the case of 
the linear synchronization distance (4) with d\ = d. However, if the parameter d\ is 
chosen smaller than d in (4), the line L intersects the line U before the jam density 
Pmax is reached (Fig. 2(b)). In this case, if the difference d — di is chosen in a proper 
way, the fluctuations in the model with linear synchronization distance (4) may be as 
simple as for the non-linear D (5), in order to lead to qualitatively the same features. 
We also did simulations (not shown in this paper), where we replaced d in (5) by 
di < d: The qualitative results remain unchanged. 

2.4- Fluctuations of acceleration and deceleration 

In order to show the power of the basic model (1), (3) [24], the remaining 
model specifications (free and safe speeds, fluctuations) will be the same as 
introduced in different Nagel-Schreckenberg CA-models in the fundamental diagram 
approach [36, 38, 39, 22, 23, 49, 50] (with a slightly more general modeling of 
fluctuations). Nevertheless, it will be shown that all features of congested patterns 
which spontaneously occur upstream of the on-ramp as well as of their evolution (when 
the flow rate to the on-ramp is changing) arc different in our CA-models. This proves 
that the basic rules (1), (3) [24] place our CA-models in the class belonging to the 
three-phase traffic theory. 
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In particular, as in the NaSch CA-models [36, 38, 39, 22, 23, 49], the accelerations 
and decelerations are stochastic. They are implemented like in [39, 22]: In a first step, 
a preliminary vehicle speed of each vehicle v n +i is 

Vn+l = V n+ l, (11) 

where v n +i is calculated based on the system of the dynamical equations (1) - (3). In 
a second step, a fluctuation arrj n (to be specified below) is added to the value v n +i 
calculated from the first step. Finally the speed v n +i at the time n + 1 is calculated 

by 

Vn+l = max(0, min(w„ + i + aTr/ n ,v n + aT,Vf rcc ,v s , n )). (12) 

This means, that the stochastic contribution arr/n may neither lead to a speed smaller 
than zero, nor to a speed larger than what the deterministic acceleration a would give, 
taking the limitations by Wf roo and v s , n into account. 
We implement the fluctuation r\ n in (12) as 

{-I \ir<p h , 
I if p b < r < p b +p a , (13) 
otherwise 

where r denotes a random number uniformly distributed between and 1. This is 
a generalization of the random deceleration (with probability pb) in NaSch cellular 
automata models [22, 23], because with probability p a also a random acceleration can 
occur. p a + pb < 1 must be fulfilled. In a different way and for a different purpose a 
random acceleration was also introduced in CA-models by Brilon and Wu [50] . As in 
a NaSch C A-modcl in the fundamental diagram approach [39] , the probability pb in 
(13) is taken as a decreasing function of the vehicle speed v n : 

where p and Po > p are constants. This corresponds to the slow-to-start rules first 
introduced by Takayasu and Takayasu [51] and later used in the NaSch CA model 
by Barlovic et al. [39]: Vehicles escape at the downstream front of a wide moving 
jam with the mean delay time t^i = r/(l — po). As in [23] this provides the jam 
propagation through free and synchronized flows with the same velocity v g of the 
downstream jam front that corresponds to a qualitative theory and to the related 
formula v g = — l/(/O m axTdei) from [6]. p m ax = l/d is the density inside the jam. 

The probability p a of the random acceleration in (13) is also taken as a decreasing 
function of the vehicle speed v n : 

P,{v n ) = {^ US) 
1 Pa2 if V n >V p . 

where w p , p a \ and p a 2 < Pal are constants. This simulates the effect that the vehicle 
moving at low speed in the dense flow tends to close up to the leading one. Indeed, 
according to (l)-(4), (12), (13), if the probability p a of the acceleration is high, the 
effect of adapting one's speed to the speed of the leading vehicle is weak: With 
probability p a the vehicle does not reduce its speed, and it may do so until it reaches 
the minimal safe gap. Note that the tendency to minimize the space gap at low speed 
can lead in particular to the 'pinch' effect in synchronized traffic flow [4, 6], i.e., to the 
self-compression of the synchronized flow at lower vehicle speed with the spontaneous 
emergence of moving jams. 
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The tendency to minimize the space gap at low speed is automatically built into 
our models, if the lower boundary L approaches the upper line U as in Fig. 2(a) and 
(b), because then the synchronization distance is no longer larger than the security 
gap at small speed. Indeed it turns out that the speed dependence (15) of p a is not 
required in these cases for a realistic modelling. Therefore we choose the probability 
p a of acceleration in (13) to be a constant in the model variants of Fig. 2. 

2.5. Cellular automata models with cruise control within three-phase traffic theory 

Nagel and Paczuski [37] proposed a variant of the NaSch CA-model where fluctuations 
are turned off for the vehicle speed v n = ff rC c- This variant has been called the NaSch 
CA-model with cruise control [1, 37]. 

In the case of such cruise control, i.e., when fluctuations are turned off for the 
vehicle speed v n = Vf ICC , simpler CA models can be used, if either the synchronization 
distance is the non- linear one, (5) (Fig. 2 (a)), or if in the synchronization distance 
(4) di < d (Fig. 2 (b)). 

In these cases, the CA-models with cruise control within three-phase traffic theory 
consist again of the formulas (l)-(3) and (5) (or (4) where d\ < d). The formula (12) 
for the incorporation of fluctuations into the final value of the speed v n +i at time 
(n + l)r is also valid, but can be simplified, if only deceleration noise like in the 
NaSch-type cellular automats [39, 22, 23] is implemented: 

if r < Ph , 

otherwise v ' 



with the probability [37, 39] 



if v n = 

Pb(v n ) = { P if < v n < tw (17) 

if V n = Vf Iee . 

Then v n +i + aTrj n < v n+1 so that (12) can also be written in the simpler form 

v n+ i = max(0, v n +i + arT] n ). (18) 



2. 6. Summary of the new models and their parameters 

In the following sections we shall discuss the congestion patterns obtained from 
simulations of four CA-models belonging to the class of three-phase traffic theory 
as introduced above. For easier reference we specify them in Table 1 with the 
abbreviations KKW-1 to KKW-4. In addition we provide two tables containing a 
list of symbols (Table 2) and typical values for the parameters (Table 3). 

A few remarks on the time and space discretization are in order. The usual 
choice of the time step r in CA-models of traffic is r = Is [2, 36]. Like in [22] we use a 
small-scale discretisation of space: The length of cells is chosen equal to 5x = 0.5 m. 
This leads to a speed discretisation in units of dv = 1.8 km/h. Hence two vehicles 
are only considered as moving with different speeds, if this difference is equal to (or 
larger) than Sv. In the model [24] with continuous changes in the vehicle speed, two 
vehicle speeds are considered as different, if their difference exceeds a much smaller 
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Table 1. Definition of four CA-models KKW-1 - KKW-4 



Dynamical part of all KKW-models 



v n+ i = max(0, min(w frcc , v B>n , u c , n )), g n = xe, n - x n - d, 

v n + ar for g n > D n - d, 

v n + ar sgn(t^, n - v n ) for g n < D n - d, 
Wf roo , d, t and a are constants. 



Vs,n — 9n/l~ 7 Vc.n — 



Stochastical part of all KKW-models 



v n +i = max(0, min(w„+i + arr] n ,v n + ar, V{ Tee , v a , n )), x n+ \ = x n + v n+ iT, 

-1 ifr<p b , 
rj n = { 1 if p b < r < p b +p a , 
otherwise. 

Specifications of synchronization distance D n and noise rj n 



for CA-model KKW-1 (cf. Fig. 1(b)): 



D -d + kv t vJv ) - { Po iiVn = r> (v ) - S Pal if Vn < Vp 

constant parameters: k, po, p, p a i, p a 2, v p . 

for CA-model KKW-2 (cf. Fig.2(a)): 



D n = d + v n r + [3vl/2a, Pb(v n ) 
constant parameters: /3, po, p, p a - 



po if v n = 

P if V n > ' 



for CA-model KKW-3 (cf. Fig.2(a)): 

PO if V n = 

D n = d + v n r + Pv n /2a, Pb{v n ) = { P ii < v n < v {ll 



if V n = Wfroo 



constant parameters: p a — 0, /3, po, p- 



for CA-model KKW-4 (cf. Fig. 2(b)): 



D n = d 1 + kv n r, 
constant parameters: d\ < d, k, po, p, p a - 



, . f Po if V n = 



value Sv = 10 -6 m/s. This is one of the reasons why fluctuations in cellular automata 
arc in general stronger than in the corresponding continuum model. 

For all CA-models within the three-phase-traffic theory investigated here we chose 
the probability p n = 0.425. This corresponds to a velocity v g — —15.5 km/h of the 
downstream front of a wide moving jam and an outflow q out — 1810 vehicles/h from 
a wide moving jam. Note that the flow rate q out refers to the case, where the vehicles 
reach the maximum vehicle speed v = Vf Tee , after they have escaped from the jam. 

For the sake of comparison with other traffic models, we choose occasionally in 
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Table 2. List of symbols 



T 


time discretization interval 


Sx 


space discretization length 


Sv — Sx/t 


speed discretisation unit 


n — n 1 9 


11U111DC1 (Jl L11I1C bLCpb 


V n 


t m ri i a \ a cnn/^n o i~ fimi^ c at~\ n unf n Anf tiiiPTii'itmn' nort 
VU111C1L. hp CCU. dL 11111C o LCp 1 i W 1 111U UL 11 ULL UelLlllg, JJdl L 


«n 


vehicle speed at time step n 




speed of the leading vehicle at time step n 


^s,n 


safe speed at time step ti 


^frcc 


maximal speed (free flow) 




vemcie posiLion ai liic Lime sicp it 




position of the leading vehicle at time step n 


a 


vehicle length 


Qn — %£,n <&n d 


gap (front to end distance) at time step n 


D 


syiiciiiomzaLioii oisiancc ai Lime SLep /£ 


a 


vehicle acceleration 





T TS~\ Kl 1 A 1 A A A /"* A 1 A T'rtf I An 

vemcic uccciciaiion 


Vn 


speed fluctuation at time step n 7 


Pa. 


piouauiiiiy oi vemcie acecici anon 


Pb 


probability of vehicle deceleration 


r 


random number uniformly distributed between and 1 


qo 


maximum flow rate in steady states 


Qmax 


maximum flow rate in free flow (v = Vf ICC ) 


Pmax = l/d 


density of traffic jam 


& 


velocity of downstream front of wide moving jam 


^out 


flow rate in the outflow from a wide moving jam 


Pmin 


density in free flow related to the flow rate <7 ou t 


to 


time at which vehicle influx from on-ramp is switched on 


(pinch) 
9lim 


limit flow rate in pinch region of general pattern 


(bottle) 
9out 


outflow from a congestion pattern at an on-ramp 



Sec. 5 model parameters deviating from the ones given in Table 3. In this case the 
parameter values are given in the related figure captions. 



3. Congested patterns on a homogeneous one-lane road 

All CA-modcls introduced here (Table 1) show qualitatively similar results on the 
homogeneous one-lane road. Therefore, only numerical results of a simulation of model 
KKW-1 (Fig. 1 (b)) will be presented in this section. 

For simulations of congested patterns on a homogeneous one-lane road, cyclic 
boundary conditions have been used. The one-lane homogeneous road has the length 
60000 cells (30 km). We checked that all qualitative results remain the same, if open 
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Table 3. Model parameters and characteristic values 



Common parameters and values for all KKW-models 


discretization units 


t = 1 s, Sx = 0.5 to, 
Sv — Sx/t — 1.8 km/h, a — 5v/t — 0.5 m/s 2 


model parameters 


Vhec = 108 km/h = 60 Sv, d — 7.5 to = 15 Sx, 
p = 0.425 


model results 


w g = —15.5 km/h, q out — 1810 vehicles//i, 
Pmin = 16.76 vehicles/fern, qo = 2880 vchicles//i 


CA model KKW-1 


general model parameters 


k = 2.55, v p = 50.4 km/h = 28 Sv, p ai = 0.2 


parameter-set I: 
model parameters 
model results 


p = 0.04, p a2 = 0.052 
9max ~ 2400 vehicles//i, ql?™^ ~ 1150 vehicles/ft- 


parameter-set II: 
model parameters 
model results 


p = 0.055, p &2 = 0.085 
g max w 2630 vehicles//i, g 1 ( i ^ nch) w 1000 vehicles/ft 


CA model KKW-2 


model parameters 
model results 


p = 0.04, p a = 0.052, 13 = 0.05 
9max ~ 2400 vehicles//i, 9ii^ 1 nch ' ) ~ 1150 vehicles/ft- 


CA model KKW-3 


model parameters 
model results 


p = 0.04, /3 = 0.05 
9max ~ 1460 vehicles//i, 9 1 ' h 3 ,' 1 nch ' ) ~ 1150 vehicles/ft- 


CA model KKW-4 


model parameters 


di = 2.5 m = 5 <5x, fc = 2.55, p = 0.04, p a = 0.052 



boundary conditions are used and the length of the road is large enough. 
3.1. Complex dynamics of synchronized flow 

The features of spatio-temporal pattern formation on a homogeneous road (i.e. a road 
without bottlenecks or on-ramps) are largely the same for the KKW-models considered 
here and the continuum model of Kerner and Klenov [24]. However, due to stronger 
fluctuations in the CA-models, the dynamics of perturbations is somewhat different. 

(i) If the initial flow rate g; n does not exceed a value <7 max < qo, 

< qin < (19) 

the fluctuations do not perturb the speed of a vehicle, which initially has maximal 
speed Ufroo, significantly (Fig. 3 (a-c)). In this case, the fluctuations lead to changes 
in the distances between vehicles, i.e. to a change in the vehicle density and hence to 
a change in the flow rate (black points F in Fig. 3 (c)). 

(ii) However, if 

<7max < <7in < <?0 (20) 
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the model fluctuations lead to an occurrence of inhomogeneous and non-stationary 
synchronized flow states where the vehicle speed is lower than «f roo (Fig. 3 (d- 
f)). According to Kerner's hypothesis about continuous spatio-temporal transitions 
between different states of synchronized flow in three-phase traffic theory [41, 42, 43] 
this behaviour corresponds to a complex motion within the 2D-region in the flow- 
density plane, where steady states exist (open circles S in Fig. 3 (f)). Similar 
continuous spatio-temporal transitions between different states of synchronized flow 
in agreement with Kerner's hypothesis have recently also been found in a different CA 
model by Fukui et al. [26]. 

These inhomogeneous synchronized flow states are the result of many independent 
local transitions at different road locations. As they cause a reduction of the initial 
maximal vehicle speed, they look similar to the F^S-transitions on the two-lane road 
which have been studied in [24]. However, in [24] the F^S-transition was in general 
caused by an external local perturbation, which led to the formation of a local region 
of synchronized flow (see Fig. 1 (c) in [24]). By contrast, in the CA models under 
consideration the local transitions were induced by the intrinsic model fluctuations in a 
wide range of densities given by (20). No external local perturbation was applied. As a 
consequence, a complex inhomogeneous spatio-temporal pattern of synchronized flow 
appears everywhere on the road instead of the local region triggered by a perturbation 
in the continuum model (compare Fig. 3 (d) in this article with Fig. 1 (c) in [24]). 

(iii) A homogeneous initial state with vehicle speed lower than V[ ICC , which in the 
absence of fluctuations would belong to the steady states within the 2D-region of the 
flow-density plane, remains a synchronized flow state for a long time. However, as in 
(ii) the evolution of these synchronized flow states shows a complex spatio-temporal 
behaviour of all traffic flow variables (Fig. 4) due to the intrinsic model fluctuations. 

3. 2. Emergence of wide moving jams 

The emergence of wide moving jams (Fig. 5) shows qualitatively the same features as 
in [24]. 

(i) In particular, as in [24], moving jams do not occur spontaneously, if the initial 
state with maximal vehicle speed v = Vf ree lies within the range of flow rates (19), 
where the maximal speed can be maintained. 

For a subset of these states, those with a density above a threshold for the wide 
moving jam formation, p m [ n (Figs. 3 (c, f) and 5 (c)), wide moving jams can be 
induced, but only by a very strong local perturbation. For the case when traffic flow 
with the maximal vehicle speed v = Vf TCC is formed downstream of a wide moving jam, 
the density p m i n is related to the flow rate in the outflow of the wide moving jam, q out . 

The velocity of the downstream front of the wide moving jam u g is a characteristic, 
i.e., unique, predictable and reproducible parameter which is a constant for given 
model parameters. This velocity together with the threshold point, (p m in, 9out), 
determines the characteristic line J for the downstream front of a wide moving jam 
(the line J in Figs. 3 (c, f), 4(c) and 5(c)). 

The strength of the perturbation needed to trigger a wide moving jam is highest 
for densities close to p m { n . Then it is not enough that a vehicle is forced to stop (the 
maximal amplitude of a perturbation) , but this stop must be maintained for some time 
(about 2-3 minutes at p m i n ) for a wide moving jam to nucleate in an initial traffic flow 
with maximal speed v — Vf Tee . Note that for initial densities which are only slightly 
higher than p m i n a wide moving jam often spontaneously disappears due to the high 
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fluctuations of the outflow from the jam in the KKW-1 model. 

(ii) As in [24], the line J determines the threshold of the wide moving jam 
excitation in synchronized flow: All densities in steady synchronized flows related 
to the line J are threshold densities with respect to the jam formation (the S^J- 
transition). At a given speed, the higher the density, the lower is the critical amplitude 
Sv c of a local perturbation for the S^J-transition: The critical amplitude Sv c for 
the S^J-transition reaches its maximum value at the threshold density. At a given 
difference between an initial density and the threshold density, the lower the initial 
speed, the lower the critical amplitude Sv c is. 

However, for synchronized flow states which lie in the vicinity of the safe speed (in 
the vicinity of the line U in Figs. 1 (b) and 2) the strong intrinsic model fluctuations 
in the KKW- models (Table 1) lead to the S^J-transition without the need of any 
external local perturbation. 

(iii) The latter result allows a simulation of the spontaneous emergence of wide 
moving jams (Fig. 5 (a-c)). In the initial state all vehicles move with the maximal speed 
v = Vfiee- For a flow rate in the interval (20), the initial vehicle speed v — Vf lee can 
not be maintained for a long time: Due to local transitions the vehicle speed decreases 
and states of synchronized flow are formed at some locations on the road (Fig. 5 (b, 
c) at x = 14 km) as already described in (ii) in Sect. 3.1. In these synchronized flow 
states, model fluctuations grow leading to the spontaneous emergence of wide moving 
jam (Figs. 5 (a) and (b, c) at x = 8 km). For the model parameters used in Fig. 5 
the spontaneous emergence of a wide moving jam can also occur for an initial state of 
synchronized flow with speed v- m < Vf rcc . 

Note that the fluctuation parameters p and p & 2 of both, random deceleration and 
random acceleration, are larger than in Fig. 3 and 4, where no spontaneous emergence 
of wide moving jams from synchronized flow states had been observed within the 
simulation time. Fast drivers are more "nervous" in Fig. 5 than in the previous 
figures. However, this is not the only reason for the emergence of wide moving jams 
in the present example. This can be seen by comparing the average values of the 
random contribution n to the speed for speeds larger than v p = 50.4 km/h: While 
in the previous Figures (n) = p a2 — P = 0.012, it is here more than twice as large: 
(rj) = 0.03. A positive value of (77) means that the drivers are biased towards stochastic 
acceleration rather than deceleration. A stronger bias implies a higher delay time in 
the vehicle deceleration. We will come back to the question, how this makes wide 
moving jams more likely, below, in Sec. 4.2(h). 

(iv) Thus, as in the model of Kerner and Klenov [24], in an initial traffic flow 
with the maximal speed v — Vf ICC model fluctuations can cause the spontaneous 
occurrence of synchronized flow, but not the spontaneous emergence of wide moving 
jams. The latter was only found in the KKW-models (Table 1), once synchronized 
flow was established. That synchronized flow states should occur first and only later 
the spontaneous wide moving jam, is a common feature of three-phase-traffic theory 
and agrees with empirical observations [3, 6]. Thus, in the KKW-models within the 
three-phase traffic theory the diagram of congested patterns on the homogeneous road 
is qualitatively different from those in other approaches [1, 2]. 
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4. Congested patterns at on-ramps 

Model of on- ramp 

In this section, for all simulations of congested patterns at an on-ramp a one-lane road 
of 100 km length (200000 cells) with open boundary conditions is used. The reference 
point x = is placed at the distance 20 km from the end of the road, so that it begins 
at the coordinate x = —80 km. The on-ramp starts at the point x = 16 km (32000 
cells) and its merging area was 0.3 km long (600 cells). 

For simulation of the on-ramp two consecutive vehicles on the main road within 
the on-ramp area are chosen randomly, their coordinates being denoted by x + > x~ . 
The entering vehicle is placed in the middle point between them at the coordinate 
x n = + x~ + l)/2] (here [ ] denotes the integer part), taking the speed of the 
leading vehicle v + [18]. In addition it was required that the distance between the two 
vehicles on the main road should exceed some value 

dx ( ™ in) = \v+ + 2d, (21) 

where A was chosen to be equal to 0.55, if not stated explicitely otherwise. 

Our numerical investigations have shown that the main qualitative features of 
congested patterns at the on-ramp and of the related diagrams do not change, if 
instead of these simple rules for a vehicle squeezing in from the on-ramp to the road 
more sophisticated lane changing rules are used. It is important, however, that on the 
one hand the model of the on-ramp allows a gradual change of the flow rate to the 
on-ramp q on from nearly zero to a relatively large (but for traffic flow relevant) value, 
and on the other hand does not introduce large additional speed fluctuations, when a 
vehicle from the on-ramp merges with the traffic flow on the main road. 

4-2. Diagrams of congested patterns at on-ramps 

A diagram of congested patterns at the on-ramp represents regions of spontaneous 
occurrence of congested patterns upstream of the on-ramp at different values of the 
initial flow rate to the on-ramp q on and the initial flow rate on the one-lane road q m 
upstream of the on-ramp. As will be discussed in more detail below, these regions in 
principle depend on how long one waits, i.e. in how far one samples rare events. In 
practice a typical waiting time much longer than one hour has little meaning, as real 
traffic is not a stationary stochastic process on large time scales. 

We found that the C A- models under consideration (Table 1) essentially show the 
same diagram of congested patterns which has been predicted for the three-phase 
traffic theory [8] and previously obtained in the continuous microscopic model [24]. 
However, for some parameter values there are interesting peculiarities in the KKW- 
models. 

(i) In Fig. 6 (a) and (b) the diagram of congested patterns at an on-ramp is shown 
for two different sets of the parameters of the KKW-1 model (Table 1) (parameter- 
set I and parameter-set II in Table 3, respectively). Although the KKW-1 model 
is considerably simpler than the model of Kcrner and Klcnov studied in [24], the 
main features of the diagram (Fig. 6 (a)) and the related congested patterns which 
spontaneously occur upstream of the on-ramp (Figs. 7, 8 and 9) are qualitatively the 
same as in [24]. 

( B} ( B) 

There are two main boundaries Fg and Sj on the diagram (Fig. 6). The limit 

( B) 

point of the boundary Fg at q on = is related to the maximum flow rate in free flow 



Cellular automata 



16 



where q ln — g max . The explanation of the limit point q in = g max is very simple: On 
the homogeneous one-lane road (i.e. without the on-ramp) synchronized flow occurs 
spontaneously, if the flow rate exceeds g max (range (20)). Therefore, synchronized 
patterns should spontaneously occur upstream of the on-ramp for q m = q max already 
for vanishing flow rate to the on-ramp q on — > 0. 

( B) 

Below and left of the boundary F s ' free flow occurs. Between the boundaries 

(B) ( B) 

Fg and Sj different synchronized flow patterns (SP) occur upstream of the on- 
ramp, without wide moving jams being observed. 

( B) 

Right of the boundary Sy> wide moving jams spontaneously emerge in 
synchronized flow which has been formed upstream of the on-ramp. The only difference 
compared to the results in [24] is that in Fig. 6 (a) there is no region where the 
moving synchronized flow pattern (MSP) occurs (about a possible occurrence of MSP 
see below). The diagram in Fig. 6 (a) is in accordance with general features of 
the diagrams of congested patterns at the on-ramp on a one-lane road which was 
postulated from qualitative considerations in [8] . 

( B} 

A few technical remarks about the determination of the boundary F s ' are in 
order: After the on-ramp has been switched on at t = t , there is a transient, before 
the congested patterns are established upstream of the on-ramp, and one needs some 
criterion to detect the transition from free to synchronized flow. The criterion we 
used is that the speed drops below some threshold value Vfs and stays low for at least 
4 min. The value of Vfs is chosen equal to 80 km/h. 

The delay time, until this criterion of the transition from free flow to synchronized 
flow is fulfilled, is marked as Tps in Fig. 10. This delay time can be rather short 
(« 1 — 2 min) for large values of q on (Fig. 10 (a)) and increases for decreasing 
q on (Fig. 10 (b)). For small values of q on the speed returns quickly to high values, 
whenever it happened to drop below Vps • Then one would not speak of a transition 
into synchronized flow any more (Fig. 10 (c)). This behaviour of Tps has been used 

( B} 

to determine the boundary F s ' . Scanning the (q on , (Zin)-planc on a grid the leftmost 
points were determined, where Tps is still smaller or equal 30 min (Fig. 10 (d)). 

( B) 

Similarly, a point on the boundary Sj ' is found as the leftmost point (<7on,<7in), 
where a wide moving jam emerges in synchronized flow within 60 min. Thus, the 

( B) (B) 

quantitative positions of the boundaries F s and Sj depend on the chosen time 
intervals. However, the qualitative forms of these boundaries are independent of the 
choice of the time intervals, provided they are high enough. 

(ii) The diagram of congested patterns at the on-ramp in Fig. 6 (b) is obtained for 
the same model KKW-1 (Table 1) as in (a), however with the fluctuation parameters 
parameter-set II in Table 3. As pointed out already in connection with Fig. 5, this 
means that the delay time of vehicle deceleration was increased compared to Fig. 6 

( B) ( B) 

(a). This has two effects: First, the region between the boundaries F s ' and Sj , 
where SP occur without spontaneous wide moving jam formation, is reduced. Second, 

( B) ( B) 

the boundaries F s ' and Sj ' merge in the limit point qi D — q max , where q on — > 0. 
This means that in the CA-model with higher delay time of vehicle deceleration wide 
moving jams can spontaneously occur in synchronized flow already at extremely low 
flow rates to the on-ramp in the vicinity of the point <fc n = q max . 

This can be explained by realizing that in the present model, an increase in the 
delay time of vehicle deceleration makes the effect of speed adjustment within the 
synchronization distance weaker. As a result, it becomes more likely that a vehicle 
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closes up to the leading one, i.e. reaches the minimal safe gap g n , so that model 
fluctuations more easily get amplified to cause a wide moving jam. 

(hi) In model KKW-2 (non-linear synchronization distance D(v)) the diagram of 
congested patterns at the on-ramp (Fig. 11 (a)) and the related congested patterns 
upstream of the on-ramp (Figs. 12, 13, 14) possess the same qualitative features as 
discussed in Fig. 6. 

There is a peculiarity of the diagram of congested patterns (Fig. 11 (b)) for model 
KKW-3 (cruise control) . In this case, there are no model fluctuations at the maximal 
speed v = Vf Icc . Thus, if the flow rate to the on-ramp g on = and the initial state 
is related to the maximal speed v = Vf ICC} synchronized flow with lower vehicle speed 
can not appear spontaneously up to the top flow rate = qo (see Fig. 11 (b) and 
Fig. 2(a)). 

However, already at an extremely low flow rate of q on w 1 — 2 vehicles/ h these 
synchronized flow states with lower speed spontaneously appear upstream of the on- 
ramp, because of a small disturbance of the initial flow at the on-ramp. Moreover, 
the flow rate q m — q max at which this effect occurs can be noticeably lower than q out 
in this case (Fig. 11 (b)). Therefore, at all flow rates q; n within the range 

<7max < <7in < <70 (22) 

SP occur spontaneously at the on-ramp for q on as small as 1 — 2 vehicles//i (Fig. 11 
(b))- 

We found that the KKW-4 model, which differs from KKW-2 only by having the 
synchronization distance which leads to Fig. 2 (b) rather than Fig. 2 (a), can show 
qualitative the same diagram of congested patterns as obtained for model KKW-2 
(Fig. 11 (a)). Likewise, if one replaces the non- linear synchronization distance in the 
cruise control model KKW-3 by formula (4) with di < d (as in KKW-4), the diagram 
of congested patterns can remain qualitative the same as Fig. 11 (b). 

To explain these results, first recall that in all these KKW-models the probability 
p a (15) was independent of the vehicle speed, in contrast to the previously considered 
KKW-1 model. The fact that this simplification of acceleration noise nevertheless 
allows us to simulate the qualitatively correct pattern formation in the three-phase 
traffic theory can be traced back to the difference in the 2D-regions of the steady 
states in the flow density plane for these CA-models: 

For low vehicle speeds the boundary L for these models is either tangential to the 
boundary U (Fig. 2(a)), or even coincides with it (Fig. 2(b)). On the other hand, the 
numerical study of the KKW-models shows that, if the vehicle speed in synchronized 
flow is very close to the safe speed, i.e. for synchronized states close to the boundary U, 
fluctuations easily lead to the emergence of a wide moving jam. As already mentioned, 
the purpose of introducing acceleration noise in our models is a simulation of the pinch 
effect, where moving jams emerge spontaneously. Both models in Fig. 2 are sufficiently 
sensitive to fluctuations at low speeds that the pinch effect can be achieved with a 
relatively small probability p a independent of the vehicle speed. 

In contrast, in the KKW-1 model the boundary L of the 2D-region of steady 
states is not close enough to the boundary U even at low vehicle speeds (Fig. 1 (b)). 
In particular, there are synchronized states with low vehicle speed, which lie below the 
line J. Therefore, the probability p a had to be chosen higher at low speeds in order 
to enhance the likelihood that a driver comes close to the boundary U. Then the 
pinch effect is also obtained in this model, in accordance with empirical observations 
in [6, 8]. 



Cellular automata 



18 



4-3. Synchronized flow patterns (SP) 

As in [24] , depending on the parameters cither the widening synchronized flow pattern 
(WSP), or the moving synchronized flow pattern (MSP), or the localized synchronized 
flow pattern (LSP) can occur in our CA-models. 

(i) As in [24] , WSP occurs at high initial flow rate on the road upstream of the 
on-ramp, g; n , and low flow rate to the on-ramp, q on , between the boundaries Fg ' and 

( B) 

Sj . The downstream front of WSP is localized at the on-ramp. The upstream front 
is continuously moving upstream, so that the width of WSP is gradually increasing in 
time (Figs. 7 (d, e), 9 (a, b), 12 (d), 14, 15 (b) and 16). 

(ii) Depending on the flow rates (fi n and q on the distribution of the vehicle speeds 
inside WSP can be related to states with nearly homogeneous speed (Figs. 7 (d) and 9 
(a, c)), or to non- homogeneous speed distributions, where sometimes non-stationary 
vehicle speed waves (propagating with different, negative and positive velocities) can 
occur (Figs. 7 (e) and 9 (b, d)). 

(hi) Note a peculiarity of the KKW-models: The upstream front in WSP 
which separates synchronized flow downstream and free flow upstream moves with 
a relatively high velocity (v^ 8 ^ w —40 km/h). This non-realistic velocity is a 
consequence of simplicity of the models presented. In spite of this, the models give a 
realistic qualitative description of congested patterns and their evolution. More correct 
values for the front velocity have been found in the microscopic model of Kcrncr and 
Klenov [24]. The other way to obtain realistic velocity of the front between free and 
synchronized flows may be the use of a strongly non-uniform free flow upstream of the 
on-ramp. In this case, there is a large spread of gaps between vehicles in free flow. As 
a result, the speed of the front is diminished due to the presence of vehicles with too 
small gaps between them. 

(iv) Recall that at very low flow rate to the on-ramp, g on , and high initial flow 
rate on the road upstream of the on-ramp, (fi n , in the diagram derived in [24] there is 
a region "MSP" , where the moving synchronized flow pattern (MSP) spontaneously 
occurs. In this case, after SP has emerged upstream of the on-ramp, this SP comes 
off the on-ramp and begins to move upstream. In some cases, a new SP emerges at 
the on-ramp; this SP comes off the on-ramp later, and so on. 

In contrast to the diagram of congested patterns in [24] , WSP can appear also at 
very low flow rate to the on-ramp <7 on (Figs. 6, 11). In other words, there is no region 
"MSP" in our diagrams, which look like the one postulated by Kerner [8] based on a 
qualitative consideration of the three-phase traffic theory approach for a one-lane road. 
However, in [8] it was also mentioned that fluctuations may cause MSP in the region, 
where WSP exists normally. Apparently for this reason, sometimes MSP appears at 

( B) 

very low flow rate to the on-ramp q on (close to the boundary F s ') in the region of 
WSP in the CA-models. 

This effect is shown for the KKW-3 cruise-control-model (Table 1) in Figs. 15 (d) 
and 17, where MSP usually occurs in the region marked "WSP & MSP" in the related 
diagram of congested patterns (Fig. 11 (b)) at lower flow rate to the on-ramp, q on . 
Within the region "WSP & MSP" (at a slightly higher flow rate q on than the one at 
which MSP occurs) a pattern which looks like a mixture of WSP and MSP can occur 
(Fig. 15 (c)): Near the on-ramp this pattern resembles MSP. However, upstream of the 
on-ramp the pattern more and more transforms into WSP. Note another peculiarity of 
the KKW-3 model: When WSP occurs in this model, the vehicle speeds and the flow 
rates in this WSP are usually related to points in the flow-density plane which are 
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in the vicinity or lie on the boundary L in the flow-density plane (Fig. 2 (a)), which 
corresponds to the synchronization distance D (see circles in Fig. 16 (b)). Presumably, 
this behaviour is not a common feature of WSP, but due to the simplicity of the model, 
(v) As in the diagram in [24], at higher flow rate q OD and lower flow rate (fr n 

( B) ( B} 

between the boundaries F s ' and Sj ' the localized synchronized flow pattern (LSP) 
occurs. The downstream front of LSP is localized at the on-ramp. However, the 
upstream front of LSP is not continuously moving upstream, so that the width of LSP 
remains spatially limited (Figs. 6 and 7 (f)). Note that the upstream front of LSP 
and therefore the width of LSP can oscillate in time. Moreover, at flow rates close to 
the boundary Fg , fluctuations can cause random appearance and disappearance of 
LSP. The boundary which separates the region of WSP (see Figs. 6 and 11) from the 
region of LSP is marked by the letter W in the diagram of congested patterns. 

4-4- General patterns (GP) 

Right of the boundaries Sj B ^ and G in Figs. 6 and 11 one finds the "general 
pattern" (GP). It is a self-maintaining congested pattern, where synchronized flow 
occurs upstream of the on-ramp, and wide moving jams spontaneously emerge in this 
synchronized flow (Figs. 7 (a, b) and 8, 12 (a, b) and 13, 15 (a)). In other words, 
in the GP wide moving jams are continuously generated somewhere upstream of the 
on-ramp. In the outflow of the wide moving jams either synchronized flow or free flow 
occurs. GP in the KKW-models within the three-phase-traffic theory have common 
features, which are very similar to those found in [24]: 

(i) If free flow occurs in the outflow of a wide moving jam, then the mean velocity 
of the downstream jam front v g and the mean flow rate in the jam outflow g ou t are 
characteristic quantities of the model. They do not depend on initial conditions and 
are the same for different wide moving jams. The mean velocity of the downstream 
front remains a characteristic parameter no matter, what the state of flow in the jam 
outflow is. 

(ii) If qi n > g Q ut, then it is obvious that the width of the wide moving jam, 
which is furthest upstream, increases monotonously (Figs. 7 (a), 8 (b), and 12 (a)). 
If in contrast, gi n < g ou t, the width of the most upstream wide moving jam decreases 
and this jam dissolves. This process of wide moving jam dissolution repeates for the 
next most upstream jam and so on (Figs. 7 (b) and and 12 (b)). Nevertheless, if 
the difference between g out and gi n is not very large, the region of wide moving jams 
is widening upstream over time (Figs. 7 (b) and and 12 (b)). Thus, GP which is 
very similar to the GP found in [24] spontaneously occurs in all C A- models (Tabic 1) 
within three-phase traffic flow theory under consideration. However, there are some 
peculiarities of the KKW-models which will be considered below. 

(iii) The first peculiarity of the KKW-models considered here is linked to the fact 
mentioned above that the upstream front in WSP which separates synchronized flow 
downstream and free flow upstream moves with a very high (non-realistic) negative 
velocity. Let us consider a case, in which the flow rate gi n is high (it corresponds 
to a point above the boundary W, Figs. 6 and 11) and the flow rate g on is related 

( B) 

to a point right of the boundary Sj in the diagram of congested patterns. In this 
case, first WSP occurs which further transforms into the GP, i.e., wide moving jams 
spontaneously emerge inside the synchronized flow of the initial WSP (Figs. 7 (a) 
and 12 (a)). However, the upstream front of this initial WSP moves considerably 
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faster upstream than the fronts of any wide moving jam. For this reason, the upstream 
front of the whole GP at any flow rate q m is determined by this upstream front of 
synchronized flow rather than by the most upstream wide moving jam. In Fig. 8 this 
upstream front of the synchronized flow is marked by the dashed line. Note that such 
a GP has been observed in empirical observations (see Fig. 20 and Sec. V.A in [8]). In 
all CA-models within the three-phase-traffic theory under consideration GP possesses 
similar non-linear features (compare Fig. 7 (a, b) with Figs. 12 (a, b), 15 (a) and Fig. 8 
with Fig. 13). 

(iv) There is also some difference in the GP formation in the continuous model [24] 
and in the KKW-models in the three-phase traffic theory. This difference concerns 
the boundary G which separates the dissolving general pattern (DGP) and the GP. 

( B) 

Recall that DGP appears right of the boundary Sj at the initial flow rate 
qin > <Zout- In this case, after a wide moving jam has been formed in synchronized flow 
upstream of the on-ramp, the mean flow rate in the jam outflow cannot exceed g ut- 
Thus, the initial condition q^ > q out is not fulfilled any more. As a result, the GP 
transforms into DGP, where one or several wide moving jams propagate upstream, 
and either free flow or one of SP occurs upstream of the on-ramp. This behaviour is 
realized also in the KKW-models (Figs. 7 (c) and 12 (c)). 

( B) 

In [24] the boundary G intersects the boundary Sj in the point <7; n = q OU f In 
the CA-models under consideration, however, the boundary G is shifted to the left in 
the diagram of congested patterns, i.e., the boundary G intersects the boundary Sj 
at some q- m > q out (Figs. 6 and 11). 

This behaviour may be explained by hysteresis effects or by the influence of high 
amplitude fluctuations (see the related remark at the end of Sec. VII. B.l in [8]). 
Indeed, in comparison with the model in [24] in the CA-models model fluctuations 
are very high. High amplitude fluctuations occur also in the outflow of a wide moving 

( B) 

jam. This may explain why for <7; n > q out right of the boundary Sj still the general 
pattern rather than DGP may occur at considerably lower flow rate to the on-ramp 
q on than in the model in [24] . 

(v) In some cases the GP like that shown in Fig. 7 (b), i.e., the GP where the 
most upstream jam is dissolved in the course of time, can occur even if the initial 
flow rate q ln is slightly higher than q out . This is linked to the fact mentioned in 
item (iii) above that the upstream front of synchronized flow in the GP propagates 
upstream faster than the one of any wide moving jam (see Fig. 7 (a, b) and Fig. 8 (b), 
where this upstream front is marked by the dashed line). Thus, upstream of the most 
upstream wide moving jam in GP a synchronized flow is formed. The flow rate in this 
synchronized flow, g^ yn \ is always lower than the initial flow rate q m . The latter flow 
rate is realized upstream of the upstream front of synchronized flow in GP (Fig. 7 (a, 
b)). Therefore, the flow rate downstream of the upstream front of synchronized flow 
in GP, i.e., the flow rate q^^ is the incoming flow rate for the most upstream wide 
moving jam rather than the initial flow rate <7; n . It can occur that q\^^ is also lower 
than q out . In this case the most upstream wide moving jam in GP will be dissolved. 

For the KKW-1 model (parameter-set I of Table 3) the maximal value of the flow 
rate q m at which GP of this type occurs is 1960 vehicles/h. For the KKW-2 model 
this flow rate q m is very close to the flow rate q out = 1810 vehicles/h. For the KKW-3 
model this flow rate is qi D = 2160 vehicles/h. 

(vi) In empirical observations of GP at on-ramps, it has recently been found [8] 
that GP possesses the following characteristic feature: If the flow rate to the on-ramp 
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q on is high enough, the average flow rate in the pinch region g(P inch ) (averaged over 
a time interval which is considerably larger than the time-distance between narrow 
moving jams emerging in the pinch region of GP) reaches the limit flow rate ql^^ ■ 
This means that the flow q(P inch ) docs not decrease below q^™ ch ^ even if the flow 
rate q on further increases. This case is called the "strong" congestion [8]. In the 
strong congestion condition, GP can not exist if q m < ql?™ ^ [45, 8]. Note that in 
the "weak" congestion condition which is realized at lower flow rates q on the flow rate 
g(pmch) changes noticeably when the flow rate q on is changing. 

As in the model [24], both the strong congestion and the weak one can be 
simulated in the KKW-models under consideration. In particular, GP can not exist 
if qi n < q^m^- Indeed, at high flow rates to the on-ramp q on the boundary Sj B ^ 
transforms into a horizontal line at g; n = g^" ^ (Fig. 6 and 11). For the KKW-1 
model (parameter-set I of Table 3) we obtain q ou t / qli^^ ~ 1.57. This value is also 
approximately in accordance with the empirical finding (see the empirical formula (4) 
in [8]). 

It should be noted that in the vicinity of this horizontal line on the boundary 
Sj , after the GP has been formed (precisely, when the related point in the flow-flow 

plane in Fig. 6 and 11 has been moved above and right of the boundary Sj ), the 
strong congestion in the pinch region occurs, and the flow rate of the vehicles, which 
may actually squeeze to the road from the on-ramp, can decrease in comparison with 
the initial flow rate q on . This effect has also occurred in the model [24] at the related 
high initial flow rates q on . Nevertheless, the GP remains the GP after this decrease in 
the real q on . This is due to a hysteresis effect in the flow rate to the on-ramp q on which 
accompanies the occurrence and the disappearance of the GP when the flow rate g n 
first increases and then decreases, correspondingly. However, the detailed investigation 
of hysteresis effects is out of the scope of this paper and it will be considered elsewhere. 



4-5. Probability of the breakdown phenomenon (the F^S transition) at the on-ramp 

Let us first recall, how the breakdown phenomenon looks like in the fundamental 
diagram approach. From a numerical analysis of a macroscopic traffic flow model 
within the fundamental diagram approach Kerner and Konhauser found in 1994 [33] 
that free flow is metastable with respect to the formation of wide moving jams (F^J 
transition), if the flow rate is equal to or higher than the outflow from a jam, q ou t- The 
critical amplitude of a local perturbation in an initially homogeneous free flow, which is 
needed for the F^J transition, decreases with increasing density: It is maximal at the 
threshold density p = p m in, below which free flow is stable. p m i n is the characteristic 
density in the outflow from a wide moving jam, i.e. when q = q ou t- The critical 
amplitude becomes zero at some critical density p = p CI > p m ; n , above which free flow 
is linearly unstable. Obviously the higher the amplitude of a random local perturbation 
the less frequent it is. Hence, the likelihood that the F— *J transition occurs in a given 
time interval should increase with density (or flow rate). The probability should be 
very small at the threshold density p m i n (at the threshold flow rate q ou t), and it should 
tend to one at the critical density p CI (at the related critical flow rate in free flow). 

In 1997 Mahnke et al [52, 53] developed a master equation approach for calculating 
the probability of the F^J transition on a homogeneous road (i.e. without bottleneck). 
Based on this approach Kuhne et al [54] confirmed that the probability of the F^J 
transition in the metastable region is increasing with the flow rate in free flow. They 
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applied this result to explain the breakdown phenomenon at a highway bottleneck. 
For a recent comprehensive discussion of the breakdown phenomenon in CA-models 
and in the Kraufi et al. model in the fundamental diagram approach see also [55, 56]. 
The theories in [33, 52, 53, 54, 55, 56] belong to the fundamental diagram approach. 

In contrast to these results, in Kerner's three-phase traffic theory [3, 42, 43] 
it is postulated that metastable states of free flow decay into synchronized flow 
(F^S transition) rather than wide moving jams (F^J transition). Moving jams 
can emerge spontaneously only in synchronized flow, i.e. after a sequence of F^S^J 
transitions [6]. In particular, even the upper limit of free flow, q max , is related to 
the F^S transition: In this limit point the probability of the F— >S transition should 
be equal to one whereas the probability of the emergence of a moving jam (F^J 
transition) should be very small. Thus, in this theory the breakdown phenomenon in 
free traffic is related to the F^S transition rather than to an emergence of moving 
jams. 

In the three-phase traffic theory, it is also postulated that the breakdown 
phenomenon at a highway bottleneck (i.e., due to an on-ramp) is due to a localized 
deterministic perturbation at the bottleneck [7]. Indeed, due to this perturbation the 
probability of the F^S transition should be considerably higher at the bottleneck than 
anywhere else. Thus, in this theory there is a spatially non-homogeneous distribution 
of the probability of the breakdown phenomenon (of the F— >S transition) on a road 
(per unit of space and time): At effective bottlenecks on the road the probability of 
the breakdown phenomenon should have maxima (see Fig. 1 (b) in [7]). This explains, 
why the breakdown phenomenon mostly occurs at highway bottlenecks. 

In the following we confirm these hypotheses of three-phase traffic theory for the 
KKW-models proposed in this paper. In particular we find that a localized permanent 
perturbation indeed occurs in the vicinity of the bottleneck (due to the on-ramp), 
and that it triggers the breakdown of initially free flow at the on-ramp rather than 
anywhere else. 

As in [24], at the boundary Fa the F— >S transition occurs at the on-ramp. Due 
to this transition the vehicle speed decreases sharply at the on-ramp. A sharp drop of 
vehicle speeds at the on-ramp (or at another bottleneck) is well-known from empirical 
observations (e.g., [57, 58]). Traffic engineers have called this effect "the breakdown 
phenomenon" in traffic flow. The F^S transition has the nature of such a breakdown 
phenomenon. 

In [24] it was also shown that the F^S transition is a first order phase transition: 
It reqires nucleation, i.e. the occurrence of a local perturbation in traffic flow whose 
amplitude exceeds some critical value. This critical amplitude is a decreasing function 
of the vehicle density in free flow (see the curve Fg in Fig. 1 (b) in [24]). The higher the 
amplitude of a random perturbation (fluctuation) the less likely it is, with a probability 
distribution that for very general reasons can be assumed to decay exponentially for 
large amplitudes. Therefore we expect that the probability that the F-^S transition 
happens within a given time interval increases exponentially with the vehicle density 
in free flow. 

Such behaviour of the probability of the breakdown phenomenon, i.e. the F^S 
transition, has indeed been observed empirically by Persaud et al. [58]. In the case 
of the on-ramp, however, there is already a local permanent non-homogeneity, which 
occurs due to the squeezing of vehicles into the main road. Corresponding to [7] this 
should explain why the F^S transition occurs at the on-ramp with a considerable 
higher probability than away from the on-ramp at the same flow rate. 



Cellular automata 



23 



This local permanent (deterministic) perturbation at on-ramp determines the 
character of the boundary F s [8]. The higher the flow rate to the on-ramp q on is, 
the higher is the amplitude of this permanent perturbation. Therefore, the higher the 
flow rate to the on-ramp q on is, the lower is the flow rate <fi n on the main road upstream 
of the on-ramp, at which the related critical amplitude occurs at the bottleneck: This 
may explain the negative slope of the curve Fg in the flow- flow plane in Fig. 6(a). 

However, a real local perturbation which leads to the F^S transition at the on- 
ramp has always also a random component, i.e., the real local perturbation should 
consist of two components: (a) a permanent perturbation, the amplitude of which is 
the higher, the higher the flow rate to the on-ramp q on is, and (b) a random component. 
The latter component should lead to the F— >S transition at the on-ramp with some 
probability also, if the flow rate upstream of the on-ramp, q m , and the flow rate to the 
on-ramp, <7 on , belong to points in the flow-flow plane in Fig. 6(a) which lie to the left of 
the boundary Fg B ^ , i.e. still in the free flow region. This probability should increase, 

if the flow rate q sum = q- m + q on approaches the boundary . If these assumptions 
are correct, then the probability of the F^S transition at the on-ramp must grow, if 
the flow rate upstream of the on-ramp <fi n increases at a given constant flow rate to 
the on-ramp q on , which can be related to the results of empirical observations [58]. 

To study the probability of the F— >S transition at the on-ramp (Fig. 18) in the 
KKW-1 model (parameter-set I of Table 3) a large number of runs of the same duration 
T has been studied for given flow rates q sunl and q on . At the beginning of each run 
there was free flow at the on-ramp. For each run it was checked, whether the F^S 
transition at the on-ramp occurred within the given time interval T or not. The result 
of these simulations is the number of realizations np where the F^S transition at the 
on-ramp had occurred in comparison with the number of all realizations Np. Then 

Pfs = np/Np (23) 

is the probability that the F— transition at the on-ramp in an initial free flow occurs 
during the time interval T at given flow rates g sum and q on . 

The flow rate <; sum was changed and the procedure with all realizations was 
repeated at the same flow rate to the on-ramp q on . The flow rate q sum at which 
the F^S transition at the on-ramp occurred in all realization is therefore related to 
the probability of the F^S transition Pps = 1- We found that lower flow rates g sum 
correspond to Pfs < 1. As expected above we found indeed an exponential increase of 
the probability Pfs as a function of the flow rate qr sum at a given flow rate to the on- 
ramp q OD . This confirms the above assumptions [7] about the nature of the breakdown 
phenomenon at the on-ramp. 

Fig. 18 shows that the nucleation rate is the higher the larger q on . This can 
be inferred from the fact that the range of flow rates downstream, Qsum^ 

over which 

the nucleation probability changes by a given amount, is much narrower for higher 
q on - The stronger permanent perturbation (higher q on ) acts like a bias which makes 
it easier to overcome the nucleation barrier. 

The F— >S transition leads to the occurrence of different synchronized flow patterns 
upstream of the on-ramp which have been considered above. 

4-6. The capacity drop 

Empirical observations show that the speed breakdown at a bottleneck is in general 
accompanied by a drop in highway capacity. If there is free rather than congested flow 
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upstream of the bottleneck, the highway capacity is usually higher. This phenomenon 
is called "the capacity drop" (for a review see [57]). 

As explained in the previous section, the breakdown phenomenon in three-phase 
traffic theory is a F^S transition at the bottleneck, so that the capacity drop is the 
difference between highway capacity in free flow and in a situation, where there is 
synchronized flow upstream and free flow downstream of the bottleneck [42, 7]. Thus 
the capacity drop is not determined by the outflow q out from a wide moving jam in 
contrast to the fundamental diagram approach. 

Obviously, in order to study the capacity drop one has to consider the outflow 
from a congested bottleneck (j^k° ttle ) , which is measured downstream of the bottleneck, 
where free flow conditions are reached. In this paper we consider the special example 
of an on-ramp as bottleneck. In [7] Kerner points out that q^ Uic ^ is not just a 
characteristic property of the type of bottleneck under consideration only. It also 
depends on the type of congested pattern which actually is formed upstream of the 
bottleneck. Hence, in the case of an on-ramp, gout**' "* is expected to vary with 
(<Zon, <7in)- Obviously, <7ou° ttle ' ) only limits the highway capacity, if it is smaller than the 
traffic demand upstream of the on-ramp, q sum = qi n + q OD} i.e. if the condition 

9out ttle) (9on, qin) < <7sum(<7on, 9in) (24) 

is fulfilled. Then the congested pattern upstream from the on-ramp simply expands, 
while the thoughput remains limited by q^° ttic \ For example, if the general pattern 
(GP) is formed at the bottleneck, an increase of q m does not influence the discharge 
flow rate q^° t ttle ^ . Instead, the width of the wide moving jam, which is most upstream 
in the general pattern, simply grows. 

Assuming that (24) is fulfilled, the capacity drop is given by 

Ag = q max - gout" le \ ( 25 ) 
where q mSLX denotes the highway capacity in free flow, i.e. the maximum flow rate in 
free flow. The minimum value, which q^° tttc ^ can take, if one considers all kinds of 
congested patterns upstream from a bottleneck, should be a characteristic quantity for 
the type of bottleneck under consideration. We denote this quantity by q^-^ 1 ^ • The 
maximum of q^ Uic ^ (denoted by (/max* 10 '') is predicted to be the maximum flow rate, 
which can be realized in synchronized flow [7], qmax^ — Hence, the capacity 

drop at a bottleneck cannot be smaller than 

Aq min = q max -q^- (26) 

This general picture of the capacity drop proposed by Kerner for the three-phase 
traffic theory has been confirmed by empirical investigations [7, 8] and in numerical 
studies of congested patterns at an on-ramp [24] , as well as in this paper. 

In particular, we found that for the KKW-1 model (parameter-set I of Table 3) 
gmax' ~ 2250 vehicles//!. This flow rate corresponds to WSP at q m = q max and low flow 
rate to the on-ramp (<7 on « 10 vehicles//! ). The flow rate q^iT^ ~ 1600 vehicles//! 
is related to LSP. Thus, the discharge flow rate q^ut^ m (25) can be changed in the 
range from 1600 to 2250 vehicles//?. As the capacity of the highway in free flow (for 
parameter-set I) is q max « 2400 vehicles//?,, the capacity drop (25) can vary in the 
range from 150 to 800 vehicles//!, depending on the type of congested pattern and the 
pattern parameters. In comparison with free flow conditions the capacity drops by 
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6.25% for WSP at low flow rate to the on-ramp up to 33.3% for LSP at a high flow 
rate to the on-ramp. 

To show the dependence of the capacity drop on the congested pattern type 
and the pattern parameters more clearly, let us consider a specific example: Initially 
the flow to the on-ramp, q on = 60 vehicles//?,, is low, and the flow rate upstream 
of the on-ramp, qi n = 2160 vehicles//;, is relatively high. As a result, WSP occurs 
(Fig. 6 (a)). In this case, the discharge flow rate g^° ttlc ' = 1950 vehicles//!, i.e., 
the capacity drop (25) is Aq = 210 vehicles//;.. Let us now increase traffic demand 
upstream on the on-ramp, i.e., the flow rate q on increases up to q on = 120 vehicles/ft 
but the flow rate q- m = 2160 vehicles//; remains the same. This causes the flow rate 
of synchronized flow in WSP to decrease more than the traffic demand (q in + q on ) 
increases. For this reason, although the flow rate q on increases, the discharge flow 
rate decreases: q^ Wle ^ = 1800 vehicles//;. This leads to a larger capacity drop (25): 
Aq = 360 vehicles//;. 

If now the flow rate q on is further increased up to q on — 240 vehicles//; (the 
flow rate q- m = 2160 vehicles//; remaining unchanged), WSP transforms into GP 
(Fig. 6 (a)). This leads to a further decrease in the discharge flow rate: q^ 1110 ^ = 
1700 vehicles//;, which is due to the pinch effect in synchronized flow of GP: The pinch 
effect strongly reduces the flow rate through the pinch region of the GP: The capacity 
drop (25) increases to Aq = 460 vehicles//;. 

If the flow rate q on is once more increased up to q on = 900 vehicles//; and the 
flow rate q- m remains the same, strong congestion occurs in the pinch region of the GP. 
This decreases the discharge flow rate to (7^° ttlc ' = 1630 vehicles/Zi. Thus the capacity 
drop (25) increases to Aq = 530 vehicles//;. In all considered cases the condition (24) 
is fulfilled. 



5. Discussion 

The numerical simulations of the CA-models within the three-phase traffic theory, 
which we have developed in this paper, show that the basic vehicle motion rules (1), 
(3) introduced by Kerner and Klenov in [24] for a microscopic three-phase traffic 
theory allow a variety of simple specifications for the synchronization distance, the 
fluctuations and for the vehicle acceleration and deceleration which lead to qualitatively 
the same diagram of congested patterns and to the same pattern features at on-ramps. 
This is due to the introduction of the synchronization distance D in the basic vehicle 
motion rules (1), (3), which leads to a 2D-region of steady states in the flow density 
plane. The robustness of the phenomena with respect to different specifications of 
the model details leads to the expectation that these phenomena occur generically 
independent of e.g. the different laws and driver behaviours in different countries. 

In the KKW-model formulation of three-phase traffic theory based on these basic 
vehicle motion rules (1), (3) specific functions for fluctuations and for the vehicle 
acceleration and deceleration can be very simple. Nevertheless, the main features of 
the diagrams of congested patterns, which these CA-models show, are qualitatively 
the same as recently found within the three-phase traffic theory and in empirical 
observations [45, 8]. 

In [24] it has been shown that these pattern features and the related diagram 
of congested patterns at on-ramps in the three-phase traffic theory are qualitatively 
different in comparison with the diagram by Helbing et al. [1, 10], which has been 
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derived for a wide class of traffic flow models within the fundamental diagram 
approach. In particular, in that diagram of congested patterns (congested states) near 
a boundary which separates trigger stop-and-go traffic (TSG) and oscillating congested 
traffic (OCT) (in the terminology of [1, 10]) a congested pattern which is a "mixture" of 
TSG, OCT and HCT (homogeneous congested traffic) should occur [1]. This pattern, 
which at first sight looks like GP, has been used in [1] for an explanation of the pinch 
effect in synchronized flow and for the jam emergence observed in [6]. However, this 
mixture pattern has no own region in the diagram of states in [10, 15, 1]: The pattern 
transforms into TSG, if (7 on decreases, or into OCT or else HCT, if q on increases. In 
our diagrams (Figs. 6 and 11 and in [24, 8]) there are no TSG, no OCT and no HCT. 
Instead, GP exists in a very large range of flow rates q on and q ln . At a given q- m GP 
in the three-phase traffic theory does not transform into another congested pattern, 
even if q on increases up to the highest possible values. Thus GP in the C A- models 
under consideration and in [24, 8] has a qualitatively different nature in comparison 
with the mixture of TSG, OCT and HCT in [1, 10]. Note that empirical observations 
of congested patterns at on-ramps [8] confirm the theoretical features of GP found 
within the three-phase traffic theory [24, 8], rather than the theoretical features of 
either TSG, or OST, or HCT, or else of the mixture of TSG, OST, and HCT within 
the fundamental diagram approach derived in [10, 15, 1]. 

In 2000, Knospe et al. [22] proposed a version of the NaSch model where in 
addition to the previous versions (e.g., [39]) drivers react at intermediate distances 
to speed changes of the next vehicle downstream, i.e., to "brake lights". The steady 
states of this model with "comfortable driving" belong to a fundamental diagram, 
i.e. the NaSch model with "comfortable driving" is a CA-model in the fundamental 
diagram approach. The NaSch model with "comfortable driving" has been applied 
in [22, 23, 47] for a description of three traffic phases: free flow, synchronized flow and 
wide moving jams. 

In order to compare this model with the ones investigated here, the congested 
patterns, which spontaneously occur at an on-ramp, and their evolution, when the flow 
rate to the on-ramp is changing were calculated for the Nagel-Schreckenberg CA-model 
with "comfortable driving" [22, 23, 47]. These new results will be discussed below. It 
will be shown that both the diagrams of congested patterns and the pattern features 
of the KKW-models in three-phase-traffic theory are qualitatively different from those 
obtained for the Nagel-Schreckenberg CA-model with "comfortable driving" , with one 
exception, concerning the wide moving jam propagation, which will be discussed first. 

5.1. Wide moving jam propagation 

The characteristic parameters of wide moving jam propagation which were found in 
empirical observations [59, 41, 4, 7] can be reproduced in many traffic flow models in 
the fundamental diagram approach, where they have first been predicted by Kerner 
and Konhauser in 1994 [33] (see also the later papers by Bando, Sugiyama et al. [60], 
by Kraufi et al. [46], by Barlovic et al. [39], and the reviews by Chowdhury et al. [2] 
and by Helbing [1]). As already mentioned, the slow-to-start rules [51, 39, 22] allow the 
wide moving jam propagation through different traffic states and bottlenecks keeping 
the characteristic velocity of the downstream jam front. This effect has recently been 
simulated in the NaSch model with "comfortable driving" [23]. This is in accordance 
with empirical observations and the wide moving jam definition made above [41, 4, 7]. 
Because the slow-to-start rules [39, 22] arc used in our CA-models as well, these CA- 
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models within the three-phase traffic theory also show the effect of the wide moving 
jam propagation through different congested patterns and bottlenecks (Fig. 19). In 
particular, if a wide moving jam which has been formed upstream of the on-ramp 
(the jam is marked as "foreign" wide moving jam in Fig. 19) then the jam propagates 
through the on-ramp and through GP (Fig. 19 (a)) and also through WSP (Fig. 19 
(b)) keeping the velocity of the downstream front of the jam. 

However, the effect of the wide moving jam propagation [41] as well as the other 
characteristic parameters of wide moving jams are apparently the only features which 
are the same in the fundamental diagram approach [33, 46, 39, 23] and in the three- 
phase traffic theory [4, 61]. All other known features of congested patterns which 
spontaneously occur upstream of the on-ramp and their evolution for the Nagcl- 
Schreckenberg CA-model with "comfortable driving" [22, 23, 47] are qualitatively 
different from those, which follow from the KKW-models within the three-phase- 
traffic theory, as will be shown in the next section. This is due to the principal 
difference between the non-linear features of congested traffic in the NaSch model 
with "comfortable driving" [22, 23] and in the KKW-models within the three-phase 
traffic theory presented in our paper. 

5. 2. Comparison with congested patterns in the Nagel-Schreckenberg cellular 
automata models with "comfortable driving" 

In Figs. 20 - 24 we compare the congested patterns and their evolution obtained in 
the KKW-models of three phase traffic theory with those for a Nagel-Schreckenberg 
type CA-model with "comfortable driving" with one lane, for which we use the rules 
and parameters presented in [22, 23, 47]. All calculations for the Nagel-Schreckenberg 
CA-model with comfortable driving arc made for two different models of the on-ramp: 
(1) The lane changing rules described in [23, 47] are applied, or (2) the model of the 
on-ramp described in Sect. 4.1 is used. In the latter case the distance between two 
consecutive vehicles on the main road, which permits a vehicle to enter from the on- 
ramp, is chosen as dx^™^ = 4d . This corresponds to the condition applied in [23, 47] 
that "an effective gap to the predecessor and a gap to the successor on the destination 
lane" is larger than or equal to the vehicle length d. It has been found that all features 
of the congested patterns and their evolution, when the flow rate to the on-ramp is 
increasing, remain qualitatively the same for both models of the on-ramp. Therefore, 
only one set of results for the model of the on-ramp described in Sect. 4.1 is shown in 
Figs. 20 - 24. This comparison shows the following results. 

(i) In the Nagel-Schreckenberg CA-model with comfortable driving [23] without 
on-ramps and other bottlenecks, i.e. on a homogeneous one-lane road, when the flow 
rate is gradually increasing, the free flow motion spontaneously transforms into a 
very complex dynamical behaviour at some critical flow rate on the road, g max : If 
a narrow moving jam emerges spontaneously, this jam dissolves within a very short 
time interval (about one- two time steps, i.e. 1-2 s), then a new narrow moving jam 
emerges which again dissolves quickly and so on at different locations and at different 
times. This behaviour resembles the oscillating congested traffic which was reported 
for other traffic flow models within the fundamental diagram approach by Lee et 
al. [11], Tomer et al. [18] and in the diagram of congested patterns at on-ramps by 
Helbing et al. [1, 10]. 

(ii) In the Nagel-Schreckenberg CA-model with comfortable driving on a one-lane 
road with an on-ramp, a complex oscillation pattern occurs upstream of the on-ramp 
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spontaneously for q in > q max , which is qualitatively the same as in (i): Narrow moving 
jams first emerge and then dissolve within a very short time interval at different 
highway locations and at different times. This pattern exists already for very small 
values of the flow rate to the on-ramp q on (Fig. 20). We will call this pattern in the 
Nagel-Schreckenberg C A- model with comfortable driving "oscillating moving jams" 
(OMJ) (Figs. 21 (d-f) and 22 (c, d)). 

The OMJ upstream of the on-ramp shows the same features as on a homogeneous 
road without on-ramp as mentioned above. The random emergence and dissolution of 
narrow moving jams on a short time scale, which is characteristic for OMJ (Figs. 21 
(figures right) and 22 (figures right)), is qualitatively different from the behaviour 
inside the widening synchronized flow pattern (WSP), which occurs spontaneously 
at the same parameters in the KKW-models within the three-phase traffic theory 
(Figs. 21 (figures left) and 22 (figures left)). 

Indeed, whereas inside WSP in the KKW-models vehicles can move with nearly 
constant vehicle speed (Figs. 21 (b) and 22 (b)), in the Nagel-Schreckenberg CA-model 
with comfortable driving inside OMJ the vehicles must randomly slow down sharply 
sometimes up to a stop and then accelerate within a short time scale, and so on 
(Figs. 21 (e) and 22 (d)). The inverse distance between vehicles shows the same high 
amplitude oscillating behaviour inside OMJ in the Nagel-Schreckenberg CA-model 
with comfortable driving (Fig. 21 (f)). Its amplitude is much higher than would be 
expected due to bare model fluctuations. This is the result of nonlinear amplification 
of fluctuations in the Nagel-Schreckenberg CA-model with comfortable driving. In 
contrast, in the KKW-models the inverse distance between vehicles shows only small 
changes in WSP, the amplitude of which is comparable to the bare model fluctuations. 
Thus, WSP in our C A- models within the three-phase traffic theory has a qualitatively 
different nature in comparison with OMJ in the Nagel-Schreckenberg CA-model with 
comfortable driving. 

(iii) If the flow rate q on is further gradually increasing first no transition to 
another congested traffic pattern occurs in the Nagel-Schreckenberg CA-model with 
comfortable driving: OMJ persists in some range of the flow rate g n (the region 
between the boundaries O and L which is marked "OMJ" in the diagram of congested 
patterns in Fig. 20). 

(iv) However, above some flow rate g n a widening region of very low mean vehicle 
speed (about v = 10 km/h) and very low mean flow rate (q^p ~ 480 vehicles/h) occurs 
spontaneously upstream of the on-ramp in the Nagel-Schreckenberg CA-model with 
comfortable driving (Fig. 23 (d-f) and 24 (c, d)). 

The downstream front of this congested pattern is pinned at the on-ramp and 
the upstream front is slowly moving upstream. Therefore, this patterns may be called 
"the widening pinned layer" (WPL for short). Inside WPL as well as in OMJ a 
very complex non-stationary behaviour of vehicles occurs. This random behaviour 
resembles the one in OMJ, however with two differences: (1) the maximal vehicle 
speed in WPL is much lower than in OMJ, and (2) the vehicles come much more 
frequently to a stop in WPL. Because of the extremely low mean vehicle speed and 
large fluctuations in WPL, we cannot discern a regular pattern in WPL. States in 
which vehicles stop can emerge and dissolve stochastically but in a correlated way 
at different locations in WPL with time scale of about 5-10 min. These correlations 
in WPL seem to propagate with a velocity faster than the propagation of the front 
between WPL and OMJ. 

(v) If the whole range of the flow rates on the road upstream of the on-ramp 
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and the flow rate <7 on to the on-ramp is studied, then the conclusion can be drawn that 
there are only two different congested patterns in the Nagel-Schreckenberg CA-model 
with comfortable driving at the on-ramp: The pinned layer (PL) which can be either 
widening (WPL) or localized (LPL) and OMJ. In addition a combination of WPL 
and OMJ is possible (the latter pattern is marked as "WPL & OMJ" in Fig. 23 (d-f) 
and 24 (c)). 

In particular, if the flow rate q\ n > qoMj (at q- m = qoMj the boundaries O and L 
intersect one another), then to the right of the boundary L in the diagram of patterns 
for the Nagel-Schreckenberg CA-model with comfortable driving WPL occurs at the 
on-ramp, and upstream of this widening pinned layer OMJ is realized (in the region 
marked as "WPL & OMJ" in Fig. 20). The case of such a spatial combination of 
WPL and OMJ is shown in Fig. 23 (d-f) and 24 (c)). 

The same combination of WPL and OMJ occurs spontaneously from an initial 
state of free flow, if the flow rate qi n is within the range gw < Qin < <Zomj- This occurs, 
when the flow rate to the on-ramp q on is increasing and becomes larger than the one 
which is related to the boundary O. 

However if the flow rate q; n is within the range qpL < qin < qw, then no OMJ 
occurs upstream of the WPL right of the boundary O (the region marked "WPL" in 
the diagram of congested patterns in Fig. 20). In the latter case free flow is realized 
upstream of the WPL. 

If the flow rate q- ln < gpL, then right of the boundary O the pinned layer occurs 
whose upstream front does not move continuously upstream. Indeed, in this case the 
flow rate <7; n upstream of the pinned layer is lower than the mean flow rate inside the 
pinned layer qpi,. This means that the localized pinned layer (LPL) rather than WPL 
occurs in the region marked "LPL" in the diagram in Fig. 20. The vehicle behaviour 
inside LPL is qualitatively the same as inside WPL. Note that OMJ and WPL, and 
also LPL, which occur in the Nagel-Schreckenberg CA-model with comfortable driving 
at the on-ramp, have not been observed in empirical observations [8] . 

These congested patterns in the Nagel-Schreckenberg CA-model with comfortable 
driving are qualitatively different from the ones which occur spontaneously for the 
same conditions in the KKW- models within three-phase traffic theory (Table 1). These 
differences are summarized below: 

(1) At the same given flow rate on the one-lane road upstream of the on-ramp, 
qi n , and at some flow rate q on , where OMJ is formed in the Nagel-Schreckenberg CA- 
model with comfortable driving, WSP occurs spontaneously upstream of the on-ramp 
in the KKW-models. It can be seen from Figs. 21 (a-c) and 22 (a, b), where the spatial 
vehicle speed distribution in the WSP is shown that these pattern characteristics are 
different from OMJ shown in Figs. 21 (d-f) and 22 (c, d). In particular, whereas OMJ 
in the Nagel-Schreckenberg CA-model with comfortable driving is characterized by a 
complex birth and decay of narrow moving jams, no moving jams were seen in WSP. 

(2) If the flow rate q on is further gradually increasing, WSP spontaneously 
transforms either into DGP or into CP in our CA-models within three-phase traffic 
theory. At high flow rate to the on-ramp q on the CP does not transform into any other 
kind of congested pattern: the GP remains GP no matter how high the flow rate q on 
upstream of the on-ramp is. In contrast, if the flow rate to the on-ramp <7 on increases 
for a given value q m > Qomj in the Nagel-Schreckenberg CA-model with comfortable 
driving, then to the right of the boundary L in Fig. 20 a widening pinned layer WPL 
develops upstream of the on-ramp. 

(3) In the Nagel-Schreckenberg CA-model with comfortable driving upstream of 
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WPL the oscillating moving jams (OMJ) occur spontaneously, if (fi n > <7omj (Figs. 20 
and Fig. 23 (d)). In the KKW-models within three-phase traffic theory, however, 
upstream of the most upstream wide moving jam in the GP a region of synchronized 
flow is realized, where no oscillations and no moving jams occur (Fig. 23 (a)). 

(4) In the pinch region of the GP in our CA-models within three-phase traffic 
theory narrow moving jams emerge. Some of these narrow moving jams grow and 
transform into wide moving jams spontaneously at the upstream front of the pinch 
region (Figs. 8 (a, b) and 7 (a)). These wide moving jams propagate further upstream 
without any limitation. Thus, in the KKW-models wide moving jams which possess 
the characteristic parameters mentioned above (Sect. 5.1) spontaneously occur in GP. 
In contrast, in the Nagel-Schreckenberg CA-model with comfortable driving narrow 
moving jams do not transform into wide moving jams. Instead, OMJ or WPL patterns 
are formed. In other words, in contrast to our CA-models within three-phase traffic 
theory, in the Nagel-Schreckenberg CA-model with comfortable driving no spontaneous 
emergence of wide moving jams which possess the characteristic parameters mentioned 
above (Sect. 5.1) occurs: In their model, such a wide moving jam can only be excited 
by an additional external perturbation of a very large amplitude (which was done 
in [23] in order to create the wide moving jam), i.e. by a perturbation which forces 
one of the vehicles to stop for several time steps. 

(5) The fact that no wide moving jams can spontaneously occur in the Nagel- 
Schreckenberg CA-model with comfortable driving can be seen from a comparison of 
the left pictures for GP with the right pictures for WPL in Fig. 23 and 24. In particular, 
there is a clear regular spatio-temporal structure of wide moving jams which alternate 
with the regions where vehicles move sometimes with nearly the maximal vehicle speed 
inside GP (left figures), whereas there is no such region inside WPL (right figures). 

Moreover, wide moving jams which have spontaneously occurred in GP (Fig. 23 
(a) and 24 (a)) propagate further without any limitation upstream. In contrast, in the 
Nagel-Schreckenberg CA-model with comfortable driving regions inside WPL, where 
vehicles come to a stop, can emerge and dissolve randomly during several minutes. 
These regions do not propagate through the upstream boundary of WPL. 

5.3. Conclusions about features of the KKW-models in the three-phase traffic theory 

Several CA-models belonging to three-phase traffic theory were proposed in this paper. 
Their simulation gives results that allow to draw the following conclusions: 

(i) The conditions (1), (3) from [24], where due to the introduction of the 
synchronization distance D a 2D-region of the steady states in the flow density 
plane appears, allow the formulation of several different sets of specific functions for 
fluctuations and for the vehicle acceleration and deceleration which lead to qualitatively 
the same diagram of congested patterns at on-ramps. 

(ii) In the KKW-model formulation, specific functions for fluctuations and for the 
vehicle acceleration and deceleration in the basic model (1), (3) can be much simpler 
than in [24]. Nevertheless, the main features of the diagrams of congested patterns 
are qualitatively the same as within the three-phase traffic theory in [8, 24]. 

(iii) The diagrams of congested patterns in the KKW-models within three-phase 
traffic theory have the following features which differ from the case considered in [24] . 
First, there is no region in the diagrams of patterns where the moving synchronized 
flow pattern (MSP) occurs exclusively. However, MSP can spontaneously randomly 
emerge in the region "WSP", where the widening synchronized flow pattern (WSP) 
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occurs. Second, when the initial flow rate on the road upstream of the on-ramp q in 
is higher than the mean flow rate in the wide moving jam outflow then the general 
pattern (GP) occurs at lower flow rates to the on-ramp in comparison with the case 
of the model in [24]. 

(iv) The congested patterns and the diagram of these patterns of the KKW- 
models explain empirical pattern features [8] and their evolution, when the flow rate 
to the on-ramp is changing. 

(v) Both the diagrams of congested patterns and the pattern features of the 
KKW-models, which belong to three-phase traffic theory, are qualitatively different 
from those derived by Helbing et al. for a wide class of traffic flow models in the 
fundamental diagram approach (for more details see [24]). 

(vi) The features of congested patterns, which occur upstream of the on-ramp in 
the KKW-models differ from those in Nagel-Schreckenberg CA-models, which belong 
to the fundamental diagram approach, including the one with comfortable driving 
(see Sect. 5.2). This is due to qualitatively different rules of vehicle motion of the 
basic model (1), (3) [24] in the three-phase traffic theory in comparison with Nagel- 
Schreckenberg CA-models [1, 2, 22, 23, 36, 37, 38, 39, 40, 49, 47]. 
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Figure 1. Hypothetical spatially homogeneous and time-independent states 
(steady states): (a) - for the initial NaSch CA-model [36]; (b) - the 2D-region 
for the steady states for the KKW-l-model (Table 1 and Table 3) is the same as 
in [24]. 
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Figure 2. 2D-regions supporting steady states in the flow-density plane for two 
of the KKW- models (Table 1 and Table 3): (a) - for KKW-2-model (non-linear 
dependence of the synchronization distance (5) on the vehicle speed) (b) - for 
KKW-4-modcl (linear synchronization distance (4) with di < d). 
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Figure 3. Traffic patterns on a homogeneous one lane road with periodic 
boundary conditions for KKW-l-model (parameter-set I of Table 3). (a, b, c) 
Perturbation of free flow at moderate flow rate (q = 2160 vchiclcs/h < <jr max ~ 
2400 vchicles/h) and (d, c, f) at high flow rate (q = 2842 vchiclcs/h > q ma x)- 
(a, d) - vehicle speed as function of time and distance (distance increases in 
downstream direction); (b, e) - vehicle speed and flow rate at the location 
x = 10 km as functions of time (one minute averaged data of virtual detectors); 
(c, f) - data in the flow-density plane which correspond to (b, e), respectively. 
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Figure 4. Synchronized flow behaviour on a homogeneous one lane road with 
periodic boundary conditions for KKW-l-model (parameter-set I of Table 3). (a) 
- vehicle speed as a function of time and distance; (b) - vehicle speed (left) and 
flow rate (right) at a fixed location (x = 10 km) as functions of time (one minute 
averages); (c) - data in the flow-density plane which correspond to (b). Initial 
free flow qi n = 1800 vchicles/h, initial speed v ln = 54 km/h (v ln = 15 m/s). 
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Figure 5. Wide moving jam formation on a homogeneous one-lane road with 
periodic boundary conditions for KKW-l-modcl with parameter-set II of Table 
3, where the probabilities of random deceleration, p, and random acceleration at 
high speed, p a 2, are larger than in Fig. 3. (a) - vehicle speed as function of time 
and distance; (b) - vehicle speed and flow rate at the locations x = 14 km and 
x = 8 km as functions of time (one minute averages); (c) - data in the flow-density 
plane which correspond to (b) at the location x = 14 km. Initial free flow rate 
q in = 2842 vehicles/h is larger than <j max = 2634 vehicles/h (cf. (19) and (20)). 
In (c) black points are related to the states of free flow with speed v close to Vf Iee 
(the points F), and circles are related to states of synchronized flow (the points 
S). 
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Figure 6. Diagrams of congested patterns at the on-ramp for the KKW-l-modcl. 
(a) - Parameter-set I of Table 3 as in Fig. 3. (b) - Parameter-set II of Table 3 as in 
Fig. 5. GP - general pattern, DGP - dissolving general pattern, WSP - widening 
synchronized flow pattern, LSP - localized synchronized flow pattern. 
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Figure 7. Congested patterns at the on-ramp belonging to Fig. 6 (a), (a) - 
General pattern (GP) at q ln > g ut, (b) - GP at qi n < q ou t, (c) - dissolving 
general pattern (DGP), (d, c) - widening synchronized flow patterns (WSP), and 
(f) - localized synchronized flow pattern (LSP). At to = 8 vain flow from the 
on-ramp is switched on. Single vehicle data are averaged over a space interval of 
40 m and a time interval of 1 min. The flow rates (9011,9111) are: (a) (500, 2300), 
(b) (740, 1740), (c) (105, 2400), (d) (90, 2300), (e) (90, 2160), and (f) (760, 1080) 
vehicles/h. 
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Figure 8. The general pattern (GP) (KKW-1 model parameters as in Fig. 7 (a)): 
(a) - vehicle speed (left) and flow rate (right), (b) - vehicle trajectories, (c) - the 
corresponding data in the flow-density plane for the location x = 15.8 km. (a, c) 
- One minute averaged data of virtual detectors whose coordinates are indicated 
in the related figures. In (c) black points are related to the states of free flow with 
the speed v close to the maximal one Vf Iee (the points F) and circles are related to 
states of synchronized flow (the points S). To show the spatio-temporal features of 
the GP clearly, only trajectories of every 6th vehicle are shown in (b). The dashed 
line in (b) shows the upstream front of the pattern which separates synchronized 
flow downstream from free flow upstream. 
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Figure 9. The widening synchronized flow pattern (WSP): (a) - The vehicle 
speed (left) and the flow rate (right), (c) - the corresponding data on the flow- 
density plane for WSP shown in Fig. 7 (d). The similar plots (b), (d) are for 
WSP shown in Fig. 7 (e). One minute averaged data of virtual detectors whose 
coordinates arc indicated in (a-d). In (c), (d) black points are related to the 
states of free flow with the speed v close to the maximal one i>f rco (the points F) 
and circles are related to states of synchronized flow (the points S). The KKW-1 
model parameters for (a, c) are the same as in Fig. 7 (d) and for (b, d) as in Fig. 7 
(e). 
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Figure 10. Speed breakdown due to the F— >S-transition at the location 
x = 15.8 km (0.2 km upstream of the begin of the on-ramp). Fixed flow rate 
q in = 2000 vchiclcs/h. Flow rate q on is: (a) 120, (b) 90, (c) 55, and (d) 70 
vehicles/h. Up-arrows mark the time to +Tps, at which according to our criterion 
synchronized flow is detected: The vehicle speed drops below the level 80 km/h 
(dashed horizontal line) and remains low for more than 4 min. Simulations of the 
KKW-l-modcl (parameter-set I of Table 3). 
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Figure 11. Diagrams of congested patterns at the on-ramp for the KKW- 
models with non-linear dependence of the synchonization distance (5) on the 
vehicle speed (Fig. 2 (a)): (a) - KKW-2-model, (b) - KKW-3-model (Tabic 1 
and Table 3). GP - general pattern, DGP - dissolving general pattern, WSP - 
widening synchronized flow pattern, MSP - moving synchronized flow pattern, 
LSP - localized synchronized flow pattern. Note that for the KKW-2-modcl (a) 
and the KKW-3-model (b) the value Vfs discussed in Sec. 4.2(i) is also chosen 
equal to 80 km/h. However for the KKW-3-model there is the exeption for the 
case, when qi n is close to the point of intersection of the curves L and F (Fig. 2 
(a)). Near this point the minimal speed in WSP can be greater than 80 km/h, and 
the threshold speed Vfs is chosen as the average between V{ ICC and this minimal 
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Figure 12. Congested patterns at the on-ramp for the KKW-2-model of Fig. 11 
(a): (a) - General pattern (GP) at <j; n > gout, (b) - GP at <j; n < q ou t, (c) 
- dissolving general pattern (DGP), (d) - widening synchronized flow pattern 
(WSP). Single vehicle data are averaged over a space interval of 40 m and a time 
interval of 1 rain, to = 8 rain. The flow rates (g n,9in) are: (a) (500, 2250), (b) 
(800, 1650), (c) (110, 2400), and (d) (70, 2300) vehicles/h. 
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Figure 13. The general pattern (GP) related to Fig. 12 (a) (KKW-2-modcl): 
(a) - vehicle speed (left) and flow rate (right), (b) - corresponding data in the 
flow-density plane at the location x = 16 km. One minute averaged data of 
virtual detectors, whose coordinates are indicated in (a, b). In (b) black points 
are related to the states of free flow with the speed v close to the maximal one 
V{ ICC (the points F) and circles are related to states of synchronized flow (the 
points S). 
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The widening synchronised flow pattern (WSP) 
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Figure 14. The widening synchronized flow pattern (WSP) related to Fig. 12 
(d) (KKW-2-model): (a) - vehicle speed (left) and flow rate (right), (b) - the 
corresponding data on the flow-density plane. One minute averaged data of virtual 
detectors, whose coordinates are indicated in (a, b). In (b) black points are related 
to the states of free flow with the speed v close to the maximal one Vf Iee (the points 
F) and circles are related to states of synchronized flow (the points S). 
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Figure 15. Congested patterns at the on-ramp for the KKW-3-model of Fig. 11 

(b) : (a) - General pattern (GP), (b) - widening synchronized flow pattern (WSP). 

(c) - widening synchronized flow pattern (WSP) at a lower value of the flow rate 
to the on-ramp than in (b), (d) - moving synchronized flow pattern (MSP). Single 
vehicle data are averaged over a space interval of 40 m and a time interval of 
1 rain, to = 8 rain. The flow rates (g n,9in) are: (a) (480, 2300), (b) (120, 2160), 
(c) (15, 2160), and (d) (5, 2040) vehicles/h. 
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Figure 16. The widening synchronized flow pattern (WSP) related to Fig. 15 
(b) (KKW-3- model): (a) - vehicle speed (left) and flow rate (right), (b) - the 
corresponding data on the flow-density plane. One minute averaged data of virtual 
detectors, whose coordinates are indicated in (a, b). 
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Figure 17. The moving synchronized flow pattern (MSP) related to Fig. 15 (d) 
(KKW-3- model): (a) - The vehicle speed (left) and the flow rate (right), (b) - 
the corresponding data on the flow-density plane. One minute averaged data of 
virtual detectors, whose coordinates arc indicated in (a, b). 
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Figure 18. Probability of breakdown phenomenon at the on-ramp for the KKW- 
1-modcl (parameter-set I of Table 3): (a) - Time dependence of the vehicle 
speed, when the F^S transition occurs at the on-ramp (the arrow marks the 
F^S transition). Data are one minute averages of a virtual detector located at 
x = 15.8 km (200m upstream of the start of the on-ramp merging area). The 
dashed line shows the speed level 80 km/h; the characteristic duration of a sharp 
decrease in the vehicle speed (this " breakdown" is marked by the arrow) is about 
1- 2 min in agreement with empirical observations. (<?on,9i n ) have the following 
values: (200, 1660) vchiclcs/h. (b, c) - The probability P PS that the F^S 
transition occurs at the on-ramp within To =30 min (curve 1) or already within 
To =15 min (curve 2), after the on-ramp inflow was switched on (at to = 8 min), 
versus the traffic demand upstream of the on-ramp, g S um = <?i n + <?on- Results arc 
shown for two different flow rates to the on-ramp, q on = 60 vehicles/h in (b) and 
q on = 200 vehiclcs/h in (c). The following criterion that a F^S transition has 
occurred is used: The vehicle speed just upstream of the on-ramp drops below 
the level 80 km/h and then remains at nearly the same low level for more than 
4 min (cf. (a)). The probabilities were obtained from N = 40 independent runs. 
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Figure 19. Simulation of the propagation of a "foreign" wide moving jam: (a) - 
Propagation through GP, (b) - Propagation through WSP. The on-ramp is located 
at x = 12 km, and the inflow is switched on at to = 7 vain. The flow rates 
(<?on,<?in) arc: (a) (480, 2300) and (b) (110, 2160). The initial location of the 
"foreign" wide moving jam is x = 18 km, the initial jam length is 0.7 km. Single 
vehicle data are averaged over a space interval of 40 m and a time interval of 
1 min. KKW-l-model (Table 1): Parameter k of the synchronization distance 
D n is k = 2.55 for (a) and k = 4 for (b). Other model parameters are given in 
Table 3, parameter-set I. 



Cellular automata 



51 



OMJ 




WPL 
LPL 



500 1000 
q on (veh./h) 



Figure 20. Diagram of congested patterns at the on-ramp for the NaSch CA- 
model with comfortable driving. OMJ - oscillating moving jams, WPL - widening 
pinned layer, LPL - localized pinned layer. The parameters of the model are taken 
from [22, 23] , in particular, the maximal speed ff roc = 27 m/s (ff roc = 97.2 km/h), 
the minimal distance d = 7.5 m, the time step t = 1 sec, the maximum flow rate 
qo = 2817 vchiclcs/h, the flow rate out from the jam q ut = 1580 vehicles/h. The 
other specific parameters of the CA-model with comfortable driving (notations 
as in [22, 23]) arc p d = 0.1, p h = 0.95, p = 0.5, h = 8, gap security = 3, the 
cell length is 1.5 m. For the simulation of the on-ramp in the CA-model with 
comfortable driving the distance dx^ ln ^ is chosen as dx^ ln ^ = Ad. The length 
of the road is 75 km, the on-ramp is at x = 65.4 km, the length of merging area 
0.6 km. As in Figs. 6 and 11, after the on-ramp has been switched on, there 
is a delay time for the congested pattern formation upstream of the on-ramp. 
However, here the boundaries O ("Oscillating") and L ("Layer") depend more 
strongly on the delay time which is chosen for awaiting the congested pattern to 
occur upstream of the on-ramp, after the on-ramp has been switched on. The 
boundary O is determined as the points (q nj9m)j where either OMJ or WPL or 
LPL occurred upstream of the on-ramp in an initial flow with maximum speed 
v = ffree within 20 min after the on-ramp has been switched on. The time interval 
for the determination of the boundary L was 60 min. The qualitative features of 
the diagram do not depend on these time intervals. The maximum flow rate in 
free flow q m ax = 1250 vehicles/h. 
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Figure 21. Comparison of the widening synchronized flow pattern (WSP) in our 
CA-model within the three-phase-traffic theory (a-c) with the oscillating moving 
jam pattern (OMJ) at the on-ramp in the NaSch CA-model with comfortable 
driving (d-f). (a, d) - the vehicle speed as function of distance and time; at 
to = 8 min the on-ramp inflow is switched on. (b, c) - the vehicle speed as 
function of distance at a given time, (c, f) - the inverse distance between vehicles 
as function of distance along the road. In (a-f ) single vehicle data are used. For (a- 
c), the KKW-l-model (Table 1) with the following parameters is used: k = 4.5, 
p = 0.04, po = 0-5) Vp = 14 m/s p a i = 0.5 and p a 2 = 0.05. The parameter 
of the merging condition (21) at the on-ramp is A = 0.4. For (d-f) the CA- 
model with comfortable driving given in [22, 23, 47] was used. In both models, 
t)f rcc = 27 m/s, d = 7.5 m, r = 1 sec, qo = 2817 vchicles/h, q ou t = 1580 vchiclcs/h 
(Table 2), the length of the road is 75 km, the on-ramp is at x = 65.4 km, the 
length of merging area 0.6 km. The other specific parameters of the CA-model 
with comfortable driving are the same as in Fig. 20. With these parameters one 
obtains the results in (a-c) for (qi n ,q n) = (1246,600) vchicles/h, in (d-f) for 
tampon) = (1246, 190) vehicles/h. 
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Figure 22. Comparison of the vehicle trajectories related to the patterns in 
Fig. 21: Widening synchronized flow pattern (WSP) in the KKW-l-model (Table 
1) within the three-phase-traffic theory (left), oscillating moving jams (OMJ) in 
the NaSch CA-model with comfortable driving within the fundamental diagram 
approach (right), (a) - vehicle trajectories (overview) of WSP shown in Fig. 21 
(a), (b) - vehicle trajectories inside of WSP. (c) - vehicle trajectories (overview) 
of OMJ upstream of the on-ramp shown in Fig. 21 (d). (d) - vehicle trajectories 
inside OMJ. Only trajectories of every 7th vehicle are shown. 
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Figure 23. Comparison of the general pattern (GP) in the KKW-l-model (Table 
1) within the three-phase-traffic theory (left) with the congested pattern at the 
on-ramp in the NaSch CA-model with comfortable driving within the fundamental 
diagram approach (right). The latter consists of a widening pinned layer (WPL) 
going over into OMJ further upstream, (a, d) - vehicle speed as function of 
distance and time; at to = 1 min the on-ramp inflow is switched on. (b, e) - 
vehicle speed as function of distance at a given time, (c, f) - inverse distance 
between vehicles as function of distance along the road. In (a-f) single vehicle 
data are used. The simulation parameters for both models are given in Fig. 20 
and Fig. 21. In both models ((ji n ,<?on) = (1964,600) vchiclcs/h. 
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Figure 24. Comparison of the vehicle trajectories related to the patterns in 
Fig. 23: GP in the KKW-l-model (Table 1) (figures left), congested pattern at the 
on-ramp in the NaSch CA- model with comfortable driving (figures right). Model 
specifications as in Fig. 20 and Fig. 21. (a) - vehicle trajectories (overview) of GP 
shown in Fig. 23 (a), (b) - vehicle trajectories inside the region of wide moving 
jams in GP. (c) - vehicle trajectories (overview) of WPL and further upstream 
OMJ corresponding to Fig. 23 (d). (d) - vehicle trajectories inside WPL. Only 
trajectories of every 7th vehicle are shown. 



